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! ABSTRACT '* ■ 
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A concise, systematic procedure is given for determining the 

Rayleigh-Schrodinger energies and wavefunctions of degenerate states to 
arbitrarily high orders even when the degeneracies of the various states 
■ are resolved in arbitrary orders. The procedure is expressed in terms of 
an iterative cycle in which the energy through the (2n+l)st order i is 

i expressed in terms of the partially determined wavefunction through the 

i 

n-th order. Both a direct and an operator derivation are given. The two 

• «•' •• » i 

approaches are equivalent and can be transcribed into each other. The 
direct approach deals with the wavefunctions (without the use of. formal 

i ■ ■ • » 

| operators) and has the advantage that it resembles the usual treatment of 
i non-degenerate perturbations and maintains close contact with the basic 
physics. In the operator approach, the wavefunctions are expressed in 
I terms of infinite order operators which are determined by the successive 
resolution of the space of the zeroth order functions. The operator 


treatment has some similarity to that of Choi (1959), but it is more 
closely related to that of either Hirschfelder (1969) or Silverstone- 
Holloway (1971). The operator expressions are useful for: double 

perturbations, expectation values of physical properties, and problems 
involving finite dimensional Hilbert space (for example, wavef unctions 
approximated by linear basis sets). The use of variational principles 
for degenerate perturbation problems is discussed. 

* *.* ■ • s. \> 

1 ' ‘ . i 
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The object of this article is to provide a practical method for 


systematically determining the n-th order Rayleigh-Schrodinger energies 
aU( j wavefunctions of perturbed degenerate states q to 

q q 

arbitrarily high orders, even when the degeneracies of the various states 
arc resolved in arbitrary orders. . 

We start with a set of solutions to £ h e zeroth order 

equation. Next the ^[°]q ar£ USed t0 help US flnd $[l]q and ^[l]q ’ 
a act of particular solutions to the zeroth and first order perturbation 

; i 

equations which correspond to the energies ' and £^ and ^^e the 

first order equation consistent. Then we use these particular solutions 
to help us find ^ [ 2 ] q * ~[2]q * and ^[’]q * a new set of "particular 


solutions which correspond to £ 


( 0 ) 


( 1 ) 


and e 


t2) 


and, make the 


q q ' q 

first and second order perturbation equations consistent. In this way we 
proceed in an iterative fashion. We call the relationship between the 
particular solutions of the first n equations to the particular solutions 
of the first' (n-1) equations the magic formula . Once it was disciydred, 
the completion of the formulation became relatively easy. In the early* 
stages of the iterative process, the are only partially determined. 

q - * 

At a later stage, after the degeneracies have been resolved, the 1 p 

• ^ *' *7n. 

become completely determined. However, at all stages of the calculations 

- - i 

the particular solutions of the first n equations suffice to determine 
the energy through the (2n+l)st order. (in accordance with Wigner's rule 
vhich states that if the error in the vavef unction is 0(X n+ ^) then the 
error in the energy is 0(X ) ). 1 ’ 



Let us explain the contents of this paper. In Section I, those 
features of Rayleigh-Schrodinger theory are presented which are common to 
both the degenerate and non-degenerate problems. The equations and concepts 
developed in Section I are used in the formulation of the degenerate treat- 
ment. In Section II, the special notation and concepts required for the 
degenerate problems are discussed.' The detailed treatment of the first 
and second order perturbations of degenerate states is given in Section 
III. The basic iterative pattern which appears, in the treatment of these 
low orders is then generalized in' Section IV to apply to the n-th_order. 

This basic cycle is then proved by induction. In Section V, 'it is shown 
how to achieve greater solutional economy by calculating the energy 

’ i 

through the (2n+l)st order after calculating the 4>’ s through the n-th 
•order. How both the basic and the (2n+l) procedures can be incorporated 
into cycles, suitable for computer programming, is shown in Section VI. 

The special simplifications which occur in orders higher than the one in 
which the degeneracy is completely resolved are given in Section ^11 ; Up 
to this point, it has been assumed that all of the perturbation equations 
are solved exactly. In Section VIII, the' use of variational principles 
and/or fixed linear basis sets to approximate the <t>'s and determine 

,* s * ^ u- 9 .^. *• 

approximate energies is discussed. Finally, Section IX is devoted to 
deriving the formal operator treatment of degenerate perturbations in a 
form such that our <J>'s . and energies, can be expressed in terms of the 
formal operators and the resolvents which have been used 

in previous papers of Hirschfelder and Silverstone and Holloway . These 
operator relations should be useful in developing interchange-type 
relations to use in the determination of physical properties. 



Those readers who are interested in the solution of the perturbation 
aquations but are not concerned with the derivations are advised ^o read 
Se ction IX , then skip to Section VI . 

Our original intention was to improve the operator treatment given in 
Reference (1). In this respect, we succeeded in deriving the operator 
fo'rmalism given in Section IX which is not inituitive and is more concise 
and easier to use than Reference (1) . In the first part of the operator 
treatment, the Hilbert space is first compacted into the ^[o]q s P ace to 
give the familiar ’'infinite order" perturbation relations. Then, by 
successive resolutions, the Hilbert space is further compacted to the space 
of the ^[n]q ' ^is point, our treatment resembles the implicit 

.... 3 ... 

formulation of Choi . In order to make further progress, v;e focus our 

.... . 

attention on the iterative cycle of starting with the relations after 

■ • *• ' i 

(n-1) resolutions; determining and ^[n]5r * an{ * then satisfying 

the consistency equation which leads to the n-th resolution. In this phase* 
our operator approach, except for a number of significant improvements, has 
a very strong resemblance (including notation) to that of Reference (1) . 


However, it became apparent that if we defined ^ [n ] J6 = ^n+l^[n]£ ( w ^ ere 
G^™^ is one of our operators) , we could transcribe our operator treatment 
into the direct formalism in which the only operators involved are the 
Hamiltonian and the resolvent. The direct treatment has the advantage of 
Having a very close resemblance. to the usual treatment of non-degenerate 
problems and it maintains a close contact with the basic physics. On the 
other hand, the formal operators may be useful in their own right since 
they can be split up iu many different ways and their Hermitean components 
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nay be moved around to form a variety of identities, such as the 
interchange-type relations which arise in connection with double 
perturbations. The formal operators may also be convenient when the 
wavefunctions are approximated by fixed linear basis sets [see Section 


VIII]. 

We believe that this paper is more complete and easier to use than 
previous formalisms. However, our treatment is not basically new since 
it represents a restructuring (by algebraic .manipulations) of the work of 
Van Vleck^, Lennar d- Jones ^ , Brillouin^, Kato^, Bloch^, Lowdin^, Primas^^, 

Morita"^, Epstein^, des Cloizeaux^, Kirtman"^, Klein‘S, Sack‘d , Chong^i 

2 3 18-22 '' . 
Silverstone , Choi , and many others . Indeed, there have been a , -• 

great many implicit treatments (where the energy is embedded in the . s ,l ■■ 

operators); and there have, been a great many explicit treatments of . 

degenerate perturbations. All of them aref really different ways of saying 

15 ' 

the same thing. Thus, Klein has made an extensive review of the explicit 

r • < ’ * • ** 

formulations and shown how they are interrelated. Then, too, des 
13 

Cloizeaux has shown how the implicit developments are related to the 
explicit. Sack^ goes even farther in showing how all of the operator 
developments are related to the Lagrangian expansion of implicitly defined 

. c • TV '> - 

variables and then he discusses the many interesting mathematical 

properties of the Lagrangian expansions. Thus, there is nothing basically 

new to be discovered in the Rayleigh-Schrodinger treatment of degenerate 

perturbations. As Einstein said of such problems, "It is simple » it. is 

just hard to do". Degenerate perturbations are very complex. However, as 
23 

Abell said, "A complexity is nothing more than a disarrangement of 
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; simplicities”. Thus, we hope that we have succeeded in rearranging 

| 

l the disarrangement. 

i 

After we had recast our development into its present direct form 

(without the use of formal operators) we discovered many points of 

I similarity with a preliminary research note given to us by Sambe and 
18 

i Roothaan . Their approach is similar in viewpoint to the Chong and 
I Larcher"^ method of "contralned secular equations”. The Sambe and 

Roothaan "wavef unctions" are really our particular solutions, the (j>'s , 
| but they lacked the relationship between the <£>’ s and the true 
; wavef unctions. They discovered the magic formula (quite independent of 
us) and their relations for the energies, including the (2n+l) 

i 

i development, are similar to ours. We hope that our paper will not 

J _ 

detract from theirs. :: 7- ; '! 


- >. • ... . 
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I. GENERAL FEATURES OF R-s‘ THEORY V7HICH APPLY TO BOTH DEGENERATE 
AND NON-DEGENERATE PERTURBATIONS. 

7 In Rayleigh-Schr6*dinger. perturbation theory we seek the solutions to 
the perturbed Schrodinger equations 


[H(X) - E z (X)]iP z (X) = 0 , 


Ci) 


where the Hamiltonian H (assumed to be Hermitean } y the wavefunction 
ip z , and the energy are assumed to be exptessible as power series 

in the perturbation parameter, 'i'i' • 


H 


H 


(0) 


n „ (n) 


+ I X n H 
n-1 


h - 4 0) + I A " h 

n=l 

. - 


* ,i. (n) 


F » r (0 > 

h 


+ I \ a c, M 


TV 


=1 


C 2 > 


The unperturbed Hamiltonian is assumed to possess a complete set 

of eigensolutions . . 


CH<°>-<°><> = -0 , '= 1 < : a< g* V 


(3) 


(0) 


where the (j>^ are a g^-fold ortho-normal set of degenerate eigen- 
functions corresponding to the energy level • The zeroth order 

wavefunction (which is defined as the limit of i|)^(X) as X 

tends to zero) is a linear combination of the • 
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In the present paper we fix our attention on a particular state q 
which arises from the perturbation of the degenerate (g^ > 1) energy 
level • However, our results also apply to non-degenerate (g^ = 1) 

cases. As a matter of fact, the following expressions for the pertur- 
bation equations, the perturbation energies, the perturbed wave fun ct ions , 
and the normalization are all equally applicable to -degenerate or non- 
degenerate problems. 

The Perturbation Equations . • 

Substituting the power series expansions (2) into the Schroclinger ' 
equation (1) for the state q yields ' * 


n 


I »" I S Ct) f (n ‘ t> - o , 

n=0 t“0 ■ • 


(4) 


where we use the notation 


5 <n) . H (n) - c (n) 
q 


(5) 


Assuming, that Eq. (4) is valid for all values of X within the radius 

of convergence of the series (2), the coefficient of each power of A 

in Eq. (4) must vanish. Thus, the n-th order Rayleigh-Schrodinger 

24 * 

perturbation equation, is 


(H (0 > - e <0) )4. (n> 

<! q 


h(t) . (n-t) 


- * 5 X 

t®! 4 


( 6 ) 


u 


In this manner the original Schrodihger Eq. (1) is resolved Into an 
infinite number of inhomogeneous perturbation, equations. 


e*v 
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In order that Eq. (6) be mathematically consistent, its right-hand 

side must be orthogonal to each of the g eigenfunctions . One of 

q qot 

these conditions is used to determine the n-th order energy; the other 

(g - 1) conditions are used in the determination of the » 

°q T q 

* * * 4 

T^ie second feature of Eq. (6) to observe is that insofar as the n-th 

,<n) 


order perturbation equation is concerned, i|> 


is the sum of a 


particular solution plus an arbitrary linear combination of the <{> 


However, such extra terms in ip 


(n) 


( 0 ) 

la 


affect the mathematical consistency 


of £he higher order perturbation equations and the determination of the 
^ * 
higher order wavef unctions . . ... .. 

_ > - ' . 

■ 7" 


The Perturbation Energies . 


The n-th order energy is obtained by multiplying Eq. (6) by ip 
24 ' ‘ q 

and integrating, 

^ <^ (0) |H (n) |/ 0) > + j l </ 0) |H (t) |^ n_t) > • (7) 


( 0 )* 


t=l 


.(n) 


! 'Here the e is expressed in terms of the wavefunction through the 

, *1 x: 

| • (n-l)st order. For n _> 3 , equivalent expressions requiring the 
4uaowledge of lower order wavefunctions can be obtained by algebraic 

' r manipulations making use of the Hermitean property of the H / and 

' •< 

I repeated use of the perturbation equations (6). Thus, for 0 k n - 1 

: c q n) “ < ^ 0) |H (n) |^ (0) >+ r n I' 1 <^ S) |fi (n *' S ~ t) |^ Ct) > (8) 

j H 4 4 s=0 t=0 4 4 

1 where the primed summations mean that the term s = t = 0 is omitted . ' 

i In accordance with the Wigner (2n+l) . rule, and e ^n+l) can be 

i . . q q 
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expressed in terras of the wavef unction through by using Eq. (8) 

with n replaced by 2n and k = n - 1 , or with n replaced by (2n+l) 
and k = n . Thus, 

e (2n) - V 0) lH (2n) |/ 0) > +- p p <i|; (s) |H (2n " s ' t) |^ (t) >, (9) 

9 9. 9 s**0 t=0 ' ** t 


and 


; (2n+l) a <i{) (0)j H (2n+l)j ii; (0) > + p p <^(s) j g (2n+l-s-t) >(1Q) 


n n 

s=0 t=»0 ^ ^ 

Although Eqs. (7)-(10) apply to degenerate, as well as to non-degenerate, 

• ; S 

perturbation problems, they are not directly useable until thd calcula- 
tions have proceeded to a sufficiently high order that the degeneracy 
is resolved and the , . . . are known. 

q q 

The Perturbed Wavef unctions . 

_ *«•-> t ^ '*’ * • 

The n-th order wavefunction can be expressed as a linear combination 
of the complete ortho-normal set of eigenfunctions of , 

- V k <n)> > <«> 


qa qa ' q. 


* k a * a 

where the prime indicates that k ^ q . Substituting' (11) into the n-th ; 
order perturbation Eq. (6), 

l' l <S* 0) - C V 0) ^ta ><<|, to ) !*q n)> - •? H Ct) »„ (n ' t) , (12) 


k a 


1=1 


.q 


so that 


„ - • - « 


^ka^^^ " lG q 0) " e k 0)] ~ 1< ^ ) | H (t) |^ n_t)> • (13) 



: - -•Thus, 


- R <0) I »«*<-*> H-Z ♦<?<♦<?’!*«>. o*> 


t=l 


a 


qa. T qa 1 q 


where R^ is the resolvent 

.(0) _ Yr r „(0) u(0) wA (0) 


- v ic-ci • 


a 


(15) 


The resolvent has the property that if F is any function 


-(H^ - e’°^)R^F 

q 


F - )F> • 


a 


qa qa 


(16) 


(0), 


Thus, if we know the function F » we can determine R F by solving 

•>5 

Eq. (16).-“ In this manner, if we know the wavefunction through the 
(n-l)st order, we can determine the part of which is orthogonal 

to the d>^ . The coefficients <d>^ must be determined from 

the normalization and the consistency conditions. 


The Normalization . 

The normalization of the' wavef unctions is not determined by the 

. - . ' ' ’• 

perturbation equations and remains arbitrary. In the- present paper we 
adopt intermediate, normalization so that 


>_ =>' 1 and |^ n ^> = 5 ^ 

q 1 q q q : n,o 


(17) 


However, it is easy to convert our results to any other type of 
normalization. A perturbed wavefunction Xq(^) corresponding to a 
different normalization scheme differs from our i/^(X) by a constant 



factor. Thus, X Q (^) (X)C (X) . Expanding the constant in powers of 

X . C (X) = X n C (n - » follows that the n-tn order wavef unction 
ft n=*0 q 

26 

corresponding to X 0 ^) is 


(n) 

X 

\ 


l r~'c 


(s),,(n-s) 


(18) 


g d 0 k 9 

For degenerate wavef unctions Eq. (18) can he expressed in the form of 
Eq. (36) , 


Solution of the Perturbation Equations . 

The Rayleigh-Schr&dinger perturbation theory is an iterative procedure 

■Tr 

for determining and E^ to any desired precision. We begin with the 

? 

solutions to the unperturbed problem, Eq, (3), and proceed recursively to 

< 

solve the perturbation eqations (6) , together with the normalization 

i 

conditions (17) , to successively higher orders of perturbation. The 
details of how this is done depends upon whether the state q under 

consideration is non- degenerate or degenerate: " 

• * » 

(1) The Nor.-Degenerate Case . 

Fcr the case that is hon-degenerate, the unperturbed 

Sehro'dinger Eq. (3) has only one independent solution' <J>^^ . Thus, 

** is known. As a result of the normalization condition (17) , 

the n-th order wavefunction (14) can be expressed in terms of the lower 
order wavefunctions " . 


-„(n) 

V q 


,( 0 ) 


The explicit determination of then requires the solution of a 

■q 

25 

partial differential equation of the Eq. (16) type. Eq. (19) together ! 


5 ii'-'y 


(t), f| (n-t) 
q 


(19) 



with the Eqs . (9) and (10) for £ and £ (2n+l) then forms the basis 

q q 

for the recursive procedure. 

(2) The Degenerate Case . 

For the case that is degenerate, is usually not known 

at the beginning of the calculation. The functions (j>^^ provide a basis 
for the so that 

q 


/ 0) = i /°><* (o v o) > . 

n L - 1 nr\ T nrv 1 ~ n 




a 


qc* qa 


( 20 ) 


Similarly, according to Eq. (14), the involves g^ unknown 

expansion coefficients • For each order n , starting with 

n = 0 , the normalization condition (17) provides one condition upon 
the expansion coefficients. The other . (g - 1) conditions are pro- 
vided by the mathematical consistency conditions for each order of the 
perturbation equations. Before discussing exactly how these equations 
are made consistent and solved, it is best that we introduce our notation 
and concepts. 


Symmetry Considerations 

In some cases, can be completely or partially determined by 

9 

group theoretic considerations. The unperturbed Hamiltonians have 

symmetry groups and G respectively. The functions d>^ are 

- qa 

bases for the irreducible representations of G^ ; whereas, 

and transform according to the irreducible representation 

of G . Thus, if is the projector of 


, then 



12b. 


If, now, commutes with as well as with H(A) , then: 

.( 0 ) 


d> - ’ = A r' are solutions to the zeroth 

qa q T qa 


(1) The set of functions ' = 

order Schrodinger equation’ and are orthogonal to the functions 

(2) All of the energy and overlap matrix elements 


-■ a-* h<°> 

Y qa q qa 


1 ( 0 ) 


vanish which link the states arising from the <j) with those which 


arise from the <j) 


( 0 ) 

qa 


(3) Thus, in determining E^ and \p , only the 


'zeroth order' functions d>^ heed to be considered and the functions 

qa 

.?( 0 ) . . 

()> ■ can be ignored. The hydrogen atom Stark effect provides a good 

example where A commutes with both and H(A) . Here 

q 

has spherical or parabolic symmetry while H(A). has only parabolic 

■resj 

( 0 ) 


symmetry. Thus A^ corresponds to a particular set of parabolic 


quantum numbers and A^(f)^“ / is one particular linear combination .of 
the hydrogen orbitals for a particular value of the principal quantum 
number. Thus, if one treats the hydrogen Stark effect in parabolic co- 
ordinates, the other degenerate functions (1-A^)4>^^ can be ignored 
and the problem becomes equivalent to a non-degenerate perturbation. 

Even if A^ does not commute with , symmetry considerations 

may be useful if there exists a subgroup which is common to both 


,( 0 ) 


and G . The irreducible representations of G are then sub- 


,27 


duced (or decomposed) into irreducible representations of its subgroups. 
The functions i and belong to irreducible representations 


,<q) 


of the common subgroup. G 


If A 1 is the projector cor- 


'1 * “ **ql 

/q) *.u a .u 


responding to D^' 1 ' , then A^ commutes with both H'” 7 and H(A) . 

Thus A = \p and the functions A play the role of the 

ql q q -ql T qa * 3 Y qa 

in the previous paragraph. Furthermore, the quantum number associated 



with A , remains unaffected by the perturbation (an example would be 

qi 

the magnetic quantum number in the Stark effect) . 

If all of the <b^ belong to the same g -dimensional irreducible . 
Iqa q 

representation of , the degeneracy is symmetry-related and not 

"accidental". If and belong to an m-dimensional irreducible 

representation of G » then the block of states containing q .will 
remain m-fold degenerate even after an infinite order of perturbation. 



.13 


NOTATION FOR DEGENERATE STATES. 


A clear, concise, and consistent notation is very important for 

treating degenerate perturbation problems where the pattern of degeneracy 

splitting cannot be predicted beforehand. The set of degenerate states 

having the zeroth order energy e ^ splits into subsets characterized 

by the first order energies e^p [where we might number 

e^P < e^P < ...]. These subsets are further split, in second order into 
ql q2 ... 

( 2 ) 

smaller subsets characterized by e After calculating the energy 

through the n-th order, the states are split into subsets corresponding 
to the energies , • e^P , eP"! , ...» . . . Thus, a 

particular state can be characterized by the "genealogy" qij...t... of 
its energy through the order in which its degeneracy is resolved. The 


degeneracy is resolved if either its subset is reduced so that it 

♦ « . i „ •- , 

contains only one state, or else, no further splittings occur in higher 
. 27 

orders . . < . - 
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In calculating the energy and wavefunctions of a set of perturbed 
degenerate states it may be necessary to use this complete genealogical 
nomenclature. However, in the development of the theory, our concern is 
with a particular state which we denote by q rather than by its full 
genealogical characterization. In addition, we use: 

[nlq labels subsets having . 

q q q ■ • • 

[n]k labels subsets having the same energy as q through (n-l)st 

Such states comprise the n-th class 

with respect to q . 

[n]£ labels subsets having the same energy as q through (n-l)st 

An) 


order but ^ 

k q 


(n) 

'£ 


order but £„ may be either equal or not equal •' e' 


To make our notation more concise, we adopt a vector notation with a 

set of functions denoted by a row vector . Thus, for example, the set 

with 1 < a < g is denoted by ijif?! . . 
qa — — & q ~[0]q^ 

In terms of our vector notation, we can discuss the resolution of the 

\ ( 0 ) 
degeneracies more precisely. In the first order, the function set ^[o]q 

is split, into the set $[jjq and the sets • Similarly, in the 

n-th order, the function set is split into the set and 

~{n-l]q * . ~[n]q 

the sets new functions $^j£ are or tho-normal linear 

combinations of the • Thus, 

~ [n-ljq- i 


m 


( 0 ) 


.(0) 


» \~S _ 1 V v / >1 \v/ \A0) ^ 

~ [n]J- -ln~l]q. ~[n-l]q ' ~[n]£ * 


( 21 ) 


and 


<A (°):- \A0) . 

^[nj£'^[n]£' 


6 £,£‘i 


( 22 ) 
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Since the ortho-normal functions span the same space as the ortho- 

normal , functions, 

~f.n-l]q 

( 0 ) . 


- [n-l]q 


,(0) *,(°> U<°) > ^ T JO) .AO) .,(0) . 

~[n]q ~ [n]q [n-l]q + M[n]k *[n]k ! ![n-l]q * ( " 3) 


(0) 


Furthermore, the projector onto the ^ basis has the resolution 

U< 0) ., >i*| 0) 1 , | - !^°> ><4,f°> I + I x+f°> I . (24) 

' ~ [n-l]q ~ [n-1] q ~[n]q -[njq 1 , L '-[njk ~[njk' 

,(0) 


And the projector onto the basis has the resolution 


><*«> | - i*; u ; i-* y y I- - czs? 

~ [OJq ~ [0]q ~ tn]q 2-[n]q‘ ^ £ 'l[s]k i[s]k' ' 


.<°> „a<0) 


n 


.(o) ;.x(o) 


If n > n* , then by the repeated use of Eq. (21), 


~ [n]& ~ [n ' ] [n ' ] q ^ ~ [n *+l) q ><: ~ [n *+l).q | * * *J ~ [n-1 3 q >< ~ [n-1] q ‘ ~ [n] & > * (26) 

Asa result of Eqs. (22) and (26), for all values of n and n 1 , 


7 *v7A 


.,(0) .,(0) . = ' ' * -I 


| Our conception is that the ^[n]q a ^ near combination of the 

: zeroth order wavefunctions of all of the states which have the energy E 


(27) 


■ through the n-th order. • Thus, the 

i 0) . ■ ■ 


is a basis .for expressing 


< 0) - feM 0) >. 


(28) 


Since zeroth order wavefunctions corresponding to. different states, are . 


orthogonal , 
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<♦<« !* (0> > 
I[n]k' q 


0 . 


(29) 


The other notational convention v/hich we adopt is to define 

> fPii > and as the zeroth order wavefunction of any of the 

r [n]q [n]k [nj£ — 

perturbed states whose energy genealogy corresponds to. [n]q, [n]k, or 

[n]£ . Thus, the functions are ii near combinations of all of the 

,( 0 ) ' 


functions ^ 


[n)k 


III. THE LOW ORDER DEGENERATE PERTURBATIONS. 


Having discussed the general features of the Rayleigh-Schrodinger 
theory, as well as the special notation and concepts required for 
degenerate problems, we are ready to proceed with the iterative process 
beginning with the zeroth order. . . . 

.t s~ ■ . . 

A. The Zeroth Order . • ■ >; •„ 

We can write the zeroth order perturbation equation either in terms 
(OV 

of the functions $[ojq w kich we knew/ 


• ■ s(0) .(0) _ 

- H J[0], - 0 _■ ■ • .. 

or in terras of the "correct zeroth order wavefunctions" ^[o]q » 


1 • 

(30) 


“"•Si, 


o , 


, (31), 


Both the ^[o]q an< ^ t * ie ^[0]q are ta ^ en to or tho-normal, 


<A (0) .;(0) > „ ,* 

[ 0 ] q [ 0 ] q" l * 


( 32 ) 
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All that we can conclude by comparing Eqs. (30) and (31) is that each of 


the 4 > 


(0) 

[0]q 


can be expressed as a linear combination of the <j> 


( 0 ) 


?[0)q ’ 


.( 0 ) 


.(0) ^(0) , „(0) 


♦iSj q - ♦[S] q <t[S q l*[ii q > 


(33) 


B." The First Order. 


.( a ) 


Each of the functions caa also he expressed by the notation 
^[1]£ w h ere > at this point in our development, the energy not 
determined. The first order perturbation Eq. (6) can be written 


5< 0 V°> + r H (1) - e (1) U (0 > = o' 

H V[in + [H £ [1]£ •[!)£ ° ’ 

But ij)!?!. is a linear combination of the » 

"IDS. ~[0]q 


(34) 


* 


( 0 ) 

U)A 




Substituting (35) into Eq. (34), 


(35) 


•s(°) A (l) + r H (l) ' (1). 1A (0) A0) . „(0) > 

H * 11 ]*, + [H " e |lJ* J *r0]q < ^t0] q l*Cl3* 


0 . 


(36) 


(1) The . ^[o] q Consistency; and ♦ 

Before solving Eq. (36), we must guarantee that it is consistent.' 

( 0 )* 

•; In order to do this, we multiply the equation by $[o]q and integrate. 
This results in the eigenvalue-eigenvector equation 


[e (1) 

l ~t0]q 


*< l > n<rf, (0) > 

- - e [ll*i 3 *[0]q**[l]*. 


0 , 


(37) 


where 


AD _ .,(0) ] P (D U (0) . 
S[0]q “ < ~ 1 0 ] q H l ^[0]q > * 


. ( 38 ) 
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Our attention shifts to the properties of the Hermitean matrix 
. It has a set of eigenvalues efiifi which are determined by 

*tO]q 

solving the secular equation ' 


(1) F (1) 

'~[0]q " E [l]£~ 


ie;~; - e r \;,,l| = 0 . 


(39) 


We designate <t ^[o]q^[l]£ > to be tbe ei £ envectors corresponding to a 

( 0 ) 

particular value of e? • They are defined by the matrix equation 


l ~[0}q 


(1) lW ,(0} 1,(0) > _ . 

[1]£~ ^[0]q'^{l]£ 0 *' 


(40) 


.<D 


Since *- G Hermitean, we can arrange so that these eigenvectors are 

orthogonal and have unit length. Furthermore, we can define the new set 
of functions as the li naar combination of the. 


♦ <»> - <*<« !+(?)> . 

~[1]£ ~f0lci ~rola‘~ril£ 


( 0 ) 


[0]q~t0]q'~[l]f. 


(41) 


The dimensionality of the t ^ ie same as the degree of degeneracy 

of the e []Lj£ • Because of Eq. (32) and the ortho-normality of the 
eigenvectors, the functions are Pitho-norraal in the sense that 




(42) 


Making use of Eq. (41), Eq. (40) can be restated 


~h (0) lH (1) U C0): > - fe Cl) - £ (1) K6 (0) U (0) > 

i[0]q |H '^riu • (e n u- c q y $ioia l hin > *. . (43) 


[i)i- 


’im 


[0]q ~ [1]£ 


Or using Eqs. (41) and (42), 

(0) IS (1) U (0) 


< $[l]£l H 


UWi ^ - fc (3-) .-(Dxx 

*^[1]£ " (e [l]£ £ q. )6 £,£’ * , 


. (44) 
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At this point we select one particular value of to call 


(i) 

an 

.a) 


; the other values of are designated as * Then we 

must consider two possible cases: 

(a) If e [i] 2 , non- degenerate, comparing Eqs. (37) and (40), it 

follows that = ^[l]£ * ' Furthermore, the eigenvector < ^[o]q^^[l]il > 

is uniquely determined by the eigenvector of Eq. (37) nr (40). Thus, 

is completely determined by Eq. (35). 

(b) If e [ j" j jj, * s degenerate, then comparing Eqs. (37) and (40) , it 

follows that must be some linear combination of the functions 

~[1]^ * -^bus, • . 5V 1 . 


.,.(0) = (0) (0) , CO) 

^Il]£ $[1]A ^£lJA*^[l3Jt * 


(45) 


where the expansion coefficients i^Jl]^ are not yat determined. 


(2) The First Order Equation; 4>f^in • s . - 

Now we are ready to return to first order perturbation equation. 
Substituting Eq. (45) into Eq. (34), 


+ [* (1) - 


On the basis of Eq. (43), Eq. (46) is consistent. Thus we can proceed to 
solve it. According to Eq. (14), the first order wavefunction can be 
expressed in the form - . . _ 


.( 1 ) _ ( 0 ) ( 1 ) ( 0 ) .( 0 ) .,( 0 ) . ( 1 ) 

*tl]t - R a *[l]t. + -E0]q < ~[0Jq ^[l]t > • 


( 47 ) 
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which, with the help of Eq. (45), we can rewrite as 


V 


(i) _ A (i) <rf> (o) i ,(o) ; (0) <A (o> i, (J d) > ,, R , 

mo ~ *[l]& ^ [1]A ~ [0]q ^[0]q |V [l]£ * (48) 


an 


where is defined to be 

~ l 1 J A/ 


<i> 


(1) 

ID* 


R (0) (1) (0) 

R H ?[!]£ * 


(49) 


From Eq. (16), it is clear that i s comp letely specified by the 

equation 


5 < 0 ) !Si + tH 


(1) .(1) , A (0) 


- e 


[1 


ys^LUJi 


(50) 


*-J7 


subject to the condition that 




o . . 


(51) 


To verify that Eq. (48) is indeed the solution of Eq. (46), it is only 

necessary to substitute (48) into (46). and remember Eq. (50). However, 

we should note that in Eq. (47) the expansion coefficients l^ri Vo >J 

~l.UJql.lJx, 

as well as the <( ^[i]q i ^ are not determined by the first order 

equation. . 


C. The Second Order. 


.( 0 ) 


Each of the functions can a i so be expressed by the notation 

^[2)£ where, at this point, the energy e [2]£ not ^ et determined. 
The second order perturbation Eq. (6) can be written 


S (( V 2) + H (1) ^ (1) + rn (2) S (2) 1t!; (0) - ft 

H ^ [2]£ + H V> [2]ji +[H “ e [ 2 ]*W [2]£ - 


: ( 52 ) 
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And since Eqs. (45) and (48) apply to any function ^[1]% including 


those which we label * 


,„(0) _ ,(0) ..(0) 1,(0)- 


[l]q ^[l]q IV t2]l 


(53) 


and 


'*$»■ - + '♦»!*» •» < 54 > 


Substituting (53) and (54) into Eq. (52), 


q(0) . (2) r5 (l),(l) + m (2) _(2) ‘(0) ]<A (0) ,.(0) > 

H ?[l]q ( ' e [23^2[13q- 1 ?[l)q^[2]^ 


+ s(l> 4 (0) <cb (0) .,(1) > , o . 

H ^[0]q $[0]q l4J [2]£ 0 v 


(55) 


First we must insure that Eq. (55) is mathematically consistent with 
respect to the functions ^[Qlq * However*; tbe function space of the 
^{Ojq spanue ^ ^y the together with the f° r values 

of k . . Thus, we divide the problem into two parts: first, to obtain 

consistency with respect to the and then consistency with respect 

to the cj>^ k 

(1) The |rijq Consistency; . ~[2]£ * 

/ O') jfc ~ - -• 

Multiplying Eq. (55) by ‘f’fQjq an<3 integrating, 

: fSui, - ‘ 0 ' w> 


_(2) _ <A (0) ]T}(1) a (1) . „(2).(0) : 

lllU "^[l)q lH ^[l]q +H ^tl]q >; ' 


( 57 ) 


where 
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Note that as a result of Eq. (43), the last term in Eq. (55) does not 

( 2 ) 

contribute to Eq. (56). As a result of Eq. (49)., *- s Kenaitean 

so that 


( 2 ) 


. — . , - una ) + h (2 V °> > 

S(l]q <H £{l)q + B ' 


(58) 


The matrix e)7n_ has the set °f eigenvalues which are 


( 2 ) 

~ [i]q 

determined by solving the secular equation 


( 2 ) 

■[ 2)1 


i_(2) _ 1 

■£[ 13 q e roio 1 l 


.( 2 ) 
[2}IV 


0 • 


(59) 


We designate <cj> r .. |4> r oio > to be the eigenvectors corresponding to s 
" 1 1 J q ~l^J x 
( 2 ) 

particular value of G [2]2, * ^ e y are defined by the matrix- equation 

(2) _ .(2) in ^(0) , A (0) 


E I[l)q " e [2U~ ]< ^tl]q^(2]2, > “ 0 * 


(60) 


( 2 ) 


Since e [1J<1 * S Hermitean, we can arrange' so that these eigenvectors are 
orthogonal and have unit length. Furthermore, we can define the new set 
of functions as the linear combination of the such that 

~ [2J& ~Ujq 


A ( 0 ) _ .kvj U ^ u ; 

~ [2]£ $[l]q ~ [ljq ~ [2] Jl 


(0) ^(0) 1.(0) 


( 61 ) 


Because of Eq. (42) and the ortho-normality of the eigenvectors, 


.,(0) |*(°) ^ = 1 X 

‘ $[2]A*$[2]A' l 6 £,r 


(62) 


Making use of Eqs. (58) and (61), Eq.. (60) becomes 


-(l).(l) , p (2) (2) v . (0) ,.(0) _ ' 

^ Jtljq (H “ € 12]^U]q i ^[2]i > ' ° 


( 63 ) 
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( 2 ) 

At this point we select one particular value of ^[ 2 ] 9 t0 ca ^ 

(2) (2) (2) 

; the other values of e ^2]il are designated as e [2]k ‘ Then we 

must consider two possible cases: 

( 2 ) 

(a) If c [ 2)£ as non- degenerate, comparing Eqs. (56) and (60), it 
follows that 4>r2]£ = * furthermore, the eigenvector I 1 * j r2]£ > 


;(o) 


is uniquely determined by Eq. (56) or (60). Thus ^[2]^ COTn pl ete ^y 

determined by Eq. (54). 

' ( 2 ) 

(b) If Gj2j£- is degenerate, then comparing Eqs. (56) and (60), 
it follows that ^[2]£ must ke some linear combination of the functions 

S[2)i • ThuSl . i, - • ’ ' ■ - 


„«>) t (0) <4 (0) u(0) 

*[2}l ?[2]£ f[2]Jt * 


(64) 


(2) The g^]k Consistenc y» ^[2] I * 

Multiplying Eq. (55) by and. integrating gived, as the result 

of Eq. (43), . . ! 4 ’ 


<(b (0 ) IhCDaW + «(?) (0) ><a (0) | (0) 

< ~ [ 1 3 k * H $[l)q + H ^l]q ^[Dq'^U]^ 


+ (e 


(1) 

UJk 


.(!>• 
q 


.( 0 ) 


" e: '^[Uk'^m* 


- 0 • 


(65) 


And, if we let 


£<-» 


[e 


( 1 ) 


e (« r 1 

e [Dk ] 


( 66 ) 


the consistency conditions (65) become 


<*(0) l,i,d) > = P (~1) »(2) <A (0) 1, (0) . 

£[l]k^[2jil C k- »[l]k ~[l]q'^[2]f ’ (67) 
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where we have defined 


a (2) 

*[!]«- 


... y ^(0) J 5 (t), .(2-t). 


Let us how express \p 


t=l 

( 1 ) 

[2J£ 


~ [ 1 3 q 


( 68 ) 


in terms of the [2)£ functions. Making 


use of Eq. (64) in the first term of Eq. (54), and using Eq. (25) to 
project the ^[Q]q functions onto the subspaces » and then using 

Eq. (67), we obtain 




( 1 ) = . .( 1 ) , ( 0 ) , ( 0 ) ■ ,( 0 ) .( 0 ) ,,.( 1 ) . 


[2]£ ' ~[l)q-£[l]q‘ v [2]£ } 

•'i 

where we have defined 


; ,(i) = r.U) + V .,(0) . nn\ 

^[2]£ { ^[l]q ^[l]q^[2]L -* „ (70) 

In Eq. (69) the functions ^[i]q ^[2)£ are com Pl ete ^y specified 

' t whereas the coefficients and <( &mt ' are as yet 


,(0) C(-1) A (2) 

'l ✓'A (0) ,,(0) 




undetermined . 


.(!> 


'.An alternate definition of the is c ^ e y satisfy the 

equation ' . ■ , r ' . . .. 


g (0) d) g (i),(0) „ 0 * 

H '?[2]£ + H t[2]Z ° *) 


(71) 


which follows from Eq. (50), subject to the orthogonality condition 


(0> ,!*«?.> ■ o , 




(72) 


which follows from Eqs. (51) and (70). A physical explanation of Eq. (72) 
is that up to this point in the development, the functions ^[2]£ are 
indistinguishable from the ^[2]£ ’ ^hus, Eq. (69) must remain valid if 
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ve replace the the corres P°nding ^[ 2 ] & functions. The 

orthogonality relation (72) makes this substituted Eq. (69) a trivial 
identity. 


A useful property of and -[2] 7 resui - ts from Eqs. (70) and 


( 68 ), 


.,(0) m (l),.(l) > . < A0) i„(2) ,.(0) . = n 

^[l]k' H i^r9ie > + ^riu-in “ 0 


L~[ 23 A . r[i]k 


12)1 .. 


(73) 


•The second order eigenvalue Eq. (60) can be restated in terms of 
and by using Eqs * ( 7 °) 311(1 (44). 


.( 2 ) _ ( ( 2 ) • ( 2 )..,( 0 ) ,.( 0 ) 
~[l)q ( (2]Jt E q > ~ [l)q ~ [2]& 


(74) 


where we used the definition (68) with 2 ■= q . , 

(2) The Second Order Equation ; ' ; ' ' * '•- 

• «...t _ - 

’ At this point we are ready to reconsider the solution of the second 

order perturbation equation. Substituting Eqs. (64) and (69) into Eq. 

(52) , 

g(0) . (2) rg (l).(l) ( (2) (2)'“ u (0) , (0) ,.(0) 

H H • *1 2]Jl (H • e [2]r*[2]£ ] $[2)£‘M2)S. 


. 5 (1),(0) .. (0) |.(1) . 

. . .^ClJ.q $ [l]q^-[2]£. 


(75) 


Using Eqs. (44), (74), and (73), it' is easy to verify that Eq. (75) is 

consistent both with respect to ancl with respect to the 1 

* 

Thus, we seek a solution. According to Eq. (14), the second order wave- 
function can be expressed in the form "... 
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(2) = r(°) ri5< 1 )^C 1 ) + (h (2) - £ (2) )h> (0} i 

[2]£ R [H M2]£ (H e [2]£^[2]Jl 

. ,(0) ..(0) . (2) 

; . + $[0]q ^[0]q‘ V t2]£ ’ 


(76) 


Using Eqs. (49), (64), and (69), Eq. (76) becomes 


.i.s**/ _ a ( 2 ) ^,( 0 ) 

^roio “ T > roio << t ) ( 


,<« 

12)1 


~ [2]£'~ [2]£ [^[2]£ > + ~ [1] q < ~ [l]q 2]£ > 


v AO) AO) . (2) 

~ [0]q $tO]q* V [2H ; 


( 2 ) 


where the $[2]l * s defined to be 


(77) 


rf,(2) r( 0) + (H^ 2 ^ - e^) 1 ' (78) 

T[2]£ U 2[2]£. ^ e [2]£ ~ [2]£ J *» Q /b) 

( 2 ) 

From Eq. (16) it is clear that 4 >j 2 ]£ is completely specified by the 
equation 


-C°) (2) g (l),<l) + m (2) -t (2) (0) 

H ^12]£ 11 2t2]£ +(H e [2]£ ) ?[2]£ 


O', (79) 


subject to the condition that 


,,(0) 1,(2) „ 

~[0]q~[2]£^ 


i ; 


(80) 


In line with the remark after Eq. (72), if the ^[2]Z functions -in Eq.. 
(77) were replaced by corresponding ^[2]£ functions, then both Eqs. 

(72) and (80) would be required in order to make this a trivial identity. 


IV. THE BASIC ITERATIVE PATTERN. 


Most of the basic iterative pattern is developed in the first dnd 
second order treatments. The remaining features appear in the third order 


J. 
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and were confirmed by considering the fourth and fifth orders. 

However, there is no need to present the third and higher orders here, 
since we prove by induction that the pattern applies to all orders. By 
pattern , we mean all of the relationships which are found in carrying out • 

i __ . 

the general n-th order iterative step. - : 

* (n) 

In the n-th order iterative step, the n-th order energies e [ n ]£ 

and the new sets of functions n n » ♦ • • r , are determined. Any 

~[n]x i[n]x, i[njx. 

zeroth order wavefuriction which corresponds to a state having the 

energy genealogy [n]Jl is a 1 linear combination of the functions ^[^3 5- * 




*<°> - U f (0 > > • 

^[n]£ . ~[n]£ -[nj^'^tnjil 


(81) 


j The first n orders of are related to the an ^ the 

^ functions determined in the previous steps by the relations 


•* 


(m) = “p 1 A (t) 

[ri)A 


~ [n-nrf-tjq “[n-m+tjq 1 / [njf, 


• .t 


(m-t). 


I - 

t « 




•C82) 


In this equation, the functions are completely specified but the 
expansion coefficients are not determined unless the -degeneracy has been 
completely resolved. • ... 

The are the particular solutions to the first n orders of 

the Rayleigh-Schrodinger perturbation equations 
. . • . m . 


5 (t) .(m-t) + ^ / (n) „(n) ^ A (0) 

~[nj£ m-,n q 


I 5 "'* 

t-0 


+ 5 > ° ' 


0<m<n ' (83) 


which satisfy the orthogonality conditions 
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“ 0 - ; 1 2 «* • 


(84) 

(85) 


The Eqs. (83)-(85) specify the to within a unitary transformation 


of the set <£ 


( 0 ) 

[n]Z • 


The wavefunctions rp 


[nJA 


are intermediately normalized in accordance 


with Eq. (17). A wavefunction Xj n j£ having a different normalization 
differs from ^[ n ]£ ' ky a constant factor, C^^^A) . Thus, substituting 


Eq. (82) into Eq. (18), the n-th order of Xj n j£ 
* • 

,(“) 

z ~ : 

m-l 


is 


Xr io = 4 > r m io < ^r°io l^r°ifl >C f°io 
*[n]& ~ [n]£ ~ [n]£' ^[n]£ [n]£ 


*» *< . . 


+ I 


(t) 


m 


l 


( 0 ) 


\* 


(s-t) >c (m-s) 


t =0 ~£ n-m + t J < l S o t ~[n-m+t]q ,v [nH [n]£ 
, (n-s) 


( 86 ) 


Since the constants C[ n j£ are determined .by the normalization scheme, 
it is easy to convert our wavefunctions into functions having a different 
normalization. 

The new sets of functions <f> r , „ are determined in terms of the 

~ [n j£ 

previously determined functions by requiring that the n-th order pertur- 
bation equation for ^[ n ]£ “^thematically consistent with respect to 
the complete set of zeroth order degenerate eigenfunctions ^[o]q * 
function space of the ^[Q]q re ?°^ ve d into the subspaces 

^[n-lJq’^tn-ljk’Jfn^Jk’' '* ’ 311(1 $[i]k ‘ 1116 consistenc y with respect 

to each subspace is considered in turn. 

The consistency requirement leads to the secular equation 

*- ' -f » * . . * - ” v * - *■ • r " * • 

for the n-th order energy €f i 0 

Injx, 
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.(n) 


I e i| 

f [n-l]q G [n]5.~ 


0 , 


(87) 


.(n) 


where e; is the Hermitean matrix 

sln-ljq 


e (n) 

~[n-l]q 


«p 


( 0 ) 


n-1 


I I + 11<n) *in-ljq > • ; < 88 > 


-- [ n- iJq V En-lJq 


(n) 

At this point we must choose which value of to call e v , the 

[n]& q 


(n) 


other values are then designated as • The eigenvectors of , 

■■■• . . . L n Jh cs L ^ T. J cj 


(0) 


are the expansion; coefficients ’ 


(n) _ (n) , (0) 

~ [n-1] q e [nJ£~ J $fn-l]q ^[n]£ r 




(89) 


The new functions are then linear combination of the , , 

~(n]£ ~[n-l]q 


,(°) 

^[n)A 


AO) U(°> > ' 

•“t n- l]q ~ [n-l]q ~ [n]f. * 


(90) 


.(0) 


■ i.V.S 


(m) 


The re q u i rement leads to the determination of the functions 


<j> ,***(« by means of ‘the magic formula, 
~ injx 


r 


j ( m ) 

J[n]l 


(*0 j. v V i(p> 


‘(-n+m-p) A (n) 


~ [n-l]q + l p | Q fln-mhplA e[n-mfp)k 




D 


‘ [n-1] q ^ [n ] Y »• - » ;; 


(91) 


! where we have defined the matrices 


A ( n >. 

zMV 


j a "Jlsir l g(t) 


; and let 


: J c - 


£(-»> . ,-(»> _ £ (u) ,-l 

C k I s , ■ [u]k J 


ff 
• / 


( 93 ) 
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We call Eq. (91) the magic formula because it is really the key to the 
general iterative procedure. 

• And, finally, the n-th order perturbation equation is solved to obtain 


4 > 


(n) 

[n)A 


T5 ( 0 ) v s(t) A (n-t) 
[n)A * 


l 


(94) 


t=l 


Eq. (94) is, of course, exactly equivalent to solving Eqs. (83) and (85) 

with m = n . Note that the functions are never required in the 

: ~ Injk 

iterative procedure for determining . 

In our derivations, we make use of the fact that when 
(n - 1) > s > p , 


(p) 


=3 <$ , e (P) _ _(p) KdT (0X 1-aW 

£[s]A' -‘Vp* [p]A r £ q } ^P^’l^p-ljq \ 


(E»- 


(95) 


The Inductive Proof of the Iterative Pattern . - : ‘ : ' ' 

‘ \ • : c . • i 

.... Our proof of the iterative scheme is inductive. We know that the 
pattern corresponding to Eqs. (81)- (95) applies through the first two 
iterative steps. Assuming that these equations apply for the first (n-1) 
steps, we prove that they apply to the n-th step. 

In treating the n-th order step, we follow the same format which we 
used for the first two steps. , In this section we focus our attention on 
a which is one of the set of j]q * Thus, in the previous 

order, the equation corresponding to Eq. (82) has n. replaced by (n-1) 
in the subscripts of the functions <p •, 


M ‘ ; °^ n - 1 • «*> 


* 
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r n \ 

Although the energy remains to be determined, the n-th order 

perturbation equation (6) can be written 


“f H (t) + <H (n) 

t“0 


(n) v . (0) = n 

" e ln]Jt )l ^[n]5, ° * 


(97) 


Substituting Eq. (96) into Eq. (97), we obtain 


H (0)7/1 (n) + r? 1 fH (n) 'p (n) U (0 > 1<d, (D) U (0) > 

H W +t i£j? ^[nrl]q +(H " e [ali ) ^a- ! lla ]< 4[n-l] B ,V [n]A* 


r C ii 


[n-l-s]g . ~[n-l-s]g |y [n]£ 


(98) 


At this point, we are ready to assure- the consistency of Eq. ^98) With 

„ „ .(0) .(0) '-.(0) ' ,(0) 

respect to the functions , f [n . 1)k , *[ n , ?]k . ••• > J U)k • 

4 * * • - > . 

(1) The Con s is t ency ; and . 

Multiplying Eq. (98) by ant * integrating, ... 


r_(n) _ (n) ,(0) ■ U<°> > „ 0 

l E[n-l]q e (n]£ i 3< ! [n-l]q 1 V [n]£ 


(99) 


where the e [^ljq ®^ ven by Eq. (88). Because of Eq. (95), the 

last line of Eq. (98) does not contribute to Eq. (99). By repeated 
manipulations involving the equations (83) [see Eq. (132) with u^n-l ). 
we can show that e [n-ijq Hermitean and it can also be written in 
the form •' ; -- > ’ - ; : ‘= : 


»[n-l]q t ^ 1 ~ [n-l]q 


~ [n-l]q ~ [n-l]q 


( 100 ) 
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As explained previously, the are the roots of the secular equation 

(87) and the corresponding eigenvectors of e j n _;jjq are the expansion 
coefficients ^[n^llq^tn]^ satisfy Eq. (89). The function 

is defined by Eq. (90). Since the eigenvectors are orthogonal 
and taken to have unit length, Eq. (84) is satisfied. 

Malting use of Eqs. (100) and (90), the conjugate complex of Eq. (88) 
becomes a new member of the Eq. (95) family. 


Also 


,(n) ' 




- ( e (n > - e (n) ')<4> (0) |d> (0) 

K [n)V, q / ; ?[n3£ ,l l[n~l]q 


(101) 


(0) i A0) _(P) 


Therefore we have proved that if Eq. (95) is valid for sjC(n-l) , then 


it is also valid for s^n . 


-• * \ 

> V*. *•' “ * 


(2) The 4> 


( 0 ) 

[n-s)k 


Consistency . - 
. ( 0 ) * 




Multiplying Eq. (98). by, ^ £ n —e Jlc anc * integrating, we obtain 


(n) 


(n) (n) w .(0) 


[A'“' ..+ )<*' w l#L> 

~[n-s]k q . [n]& ~ [n-sjk 1 ~ [n-l]q £[n~l]q 1 v [n]£ 


! , !i! p ?-> ■= ' d . 

• p Zi ~[n-sjk ~[n-l-p]q r [n]2 - 

However, by repeated use of Eq. "(90)\, 


( 102 ) 


^[n-ljq ~[nr-s]q < ~In-s]q ^~[n-s+l]q >< ~[n-s+l]q ^ * * * ^[n-ljq”* ’ ( 103 > 



• J-. 
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so that according to Eq. (84) 


<* 


( 0 ) 


I > 


~ [n-ajk 1 ~ [n-1] q 


0 . ,‘ 


Also, from Eq. (95), , 


A Ort-p ) 

~ [n-s 3k 


0 , p > s 


(n-s)- 


. ~ [n-s]k 

Thus, since by Eq. (90), 


(ef n - s > - e (n - s) )<4.f 0) ,J^ 0) > • 

[n-s]k q ~ [n-s ]k ~ [n-s-l)q ‘ 


‘.j=, : 


<tfO) . |^f S ?o> - 1 4> > <4>r 0) |t^ S ? 0 > 

~in**s]k' ~[n-s]k' ~ [n-s-l]q ~[n-s-l]q' ^[n]£. 


» 

■ ■) 


i . 


;«>) 


. • , - 


Eq. (102), the 4>j n consistency condition, becomes 


s-1 


l I > _ “y^pCs-n) (n-p) ' , (0) | . Cp> > . 

~[n-s]k‘ ^[n]X. =0^^ ~[n-s]k ~[n-l-pjq'^[n]£ 


In deriving the magic formula, these consistency conditions 
used in the following connection. Since,- 


Ao) o .(0) 'AO) |,C0) > 

~ [n-8-l]q ^[n-s]q tEn-s]q'?[n-s-l]q 


v. 


. v j,(°) AO) | oiu; 

r ~[n-s)k -[n-s]k ~[n-s-l]q 


( 0 ) 


(104) 

(105) 

(106) 

(107) ‘ 

V- V* 
i, to t • 

(108) 
are 

(109) 
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it follows that 


< f[n-s-l]q^[n]Jl > C ^[n-s-ljq ^fa-s]q > ^[n-slq^tn}? 


k ~.tn r s]k‘ [n]£* 


IaW 


(o) 


( 110 ) 


/o) 

~ [n-s]k 


Thus, using the (j) r '' / ■>, consistency condition. 


- 1 n I S T 0 > “ <<3^°^ 1 1 4> r J ^>r S T <7 > 

~[n-s-l]q ^inJJi ~ [n-s-l]q ~ [n-sjq ~ [n-s]q 1 ^£n]£ 


T <A (0) u (0) > S y 1 c(s-n) (n-p) (0) |.(p) 

£ ^[n-s-ljq^U-sJ^ p £ 0 °k £[n~s]k ^jn-l-pjq 1 V [n]il 


(111) 


\ ' 

(3) The Proof of thh Magic Formula;*. . 

-Injx, 

, We seek to show that i£ the magic formula is valid for the deter- 

(m * ) 

raination of all of the d> r for n’<n and m'<n' , then it is also 

I[n J* 

valid fjor the determination of the ^[n]S, w ^- t ^ m<n • We do this by 

comparing the two expressions for > Eqs. (96) and (82). Id 

/ in * J • *’ ' 

Eq. (96), all of the wil ich occur correspond to n'<n atii3 

therefore we can assume that they are known at this stage of the .develop 0 

/ i.\ 

ment. In Eq. (82), except for the ' ip r 1( , the only unknown function at 

Injx 

this point is * B y starting out with Eq. (96) and making ' 

repeated use of Eq. (109), we obtain a new equation which agrees with 
Eq. (82) provided that defined in accordance with the magic 

formula (91) . 


! 
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In order to make cur derivation less cumbersome, let us define the 
following functions: 


E t < P ) ’ Jtn^rf-t-Uq + 


PI 1 


u=0 


V rf (u) 

“ ~ [n-is+ulk^k 


_ (U2) 


V*> . 


and 


G { . - 


(t) 


M 


- j, \ i , h (m t) 

*X I[n-nH-t]q X[n-ia-rt3q ■ [n] Z * . . > 


we note that, according to ; the magic formula (91) , 


<113) 

(114) 


~ [n-mrftjf. r 


( 0 ) 


.(0) 




K " » 


(115) 
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Now we are ready to proceed to revise Eq. (118). V7e start by 

applying Eq. (116) Onto the term Fq( 0) • The terms in t]ie Eq. (116) 

sum provide a contribution' to each of the remaining F (0) so that 

m 

the net result of applying Eq. (116) onto Fq( 0) is that Eq. (118) 
becomes 




(m) 


m 


= G, 


|nj*. - "0 + l f t (1) •. 
t=l 


(119) 


Applying Eq. (116) upon F^(l) ; then upon F^(2) ; ... ; and finally 


on F (s) changes Eq. (119) into 
6 


* 


(m) 


s . m 

r io “ T G + l F (6+1) . 

tn)1 A v t .U - 


■ t 


• ! ( 120 ) 


When s=a-l , Eq. (120) becomes 


• -.t 


f V m-1 

io " 1 G + F (m) 

t t Q t m w 


(121) 


Comparing Eqs. (117) and (121), 


: *:v- 




( 122 ) 


Using Eq. (81) [which we have already verified], it follows that 


iw* ’ ’ 


*■(123) 


where f^On) is-. given by Eq. (112). Eq. (123) is indeed the magic 
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formula! Thus, we have completed the proof. 

(4) The_ m-th Order Perturbation Equations . 

- According to Eq. (83), the functions ^[^-ljq sat isfy the per- 
turbation equations 


m 


I 5 

t=0 


5<t) 


(m-t) 

• tn-l]q 


0<m<n-l ... 
N ** 


(124) 


Because o£ Eqs. (SI) [which we just proved] and Eq. (124), the 4> 


a («) 


[n]l 


satisfy the first (n-l)st order perturbation equations (83) , And 


r, 


making use of Eq. (91) , Eq. (89), and the orthogonality conditions (85) 
■ fm) 

for with n'<n , it follows that Eq. (85) is satisfied by the 

, (m) : - 

~[n)£ \ ‘ ' r - ■ 

» ' 

(5) , The n-th Order Perturbation Equation and * 

At this point we are ready to solve the n-th order perturbation 
equation. Substituting' the expressions Eq. (82) for with 

m<n into Eq. (97), 


«(0) . (n) 


n 


"'"♦Sii + U 


5 (t) A (n-t) 

W 


= (n) -e{ n > )4 0 h]<*;°ll^°l> 

q [n)JL~[n)& J ~[n] ^inJJl 


n-1 n-e 


+ i i r/ojn-s-o w) | (s> 

s “l ~ [n-s]q ~.[n-s]q l ^[n]£ 


(125) 


Thus, according to Eq. (14), 


, u (*0 


n 


o _(0) v ; -(t) . (n-t) . (0) | , (0) . 

J [nlt E . Jj H . 


n-1 n-s 


t i i 

s=l t“l 


C°)s(t) A (n-s-t) <({) (0) 


[n^-slq ~ [n-s]q 


'♦Si x + • 


(126) 


. J . * 
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Thus, if we use the expressions given by Eq. (94) for <j>£” S \ and 

~ tu a j q 

define in accordance with Eq. (94), then is given by 

Eq. (82). Substituting this expression for into Eq.. (125), it 

follows that the satisfy the n-th order equation (83),.- Since . 

ariy function -.of the form R^(...) is orthogonal to the ^[o]q » it: . 
follows from Eq. ( 94 ) that satisfies the orthogonality condition 

Eq. (85). * 


(6) Summary . ... •' . •.•••. 

» - - - / 

Thus we have shown that if Eq. (81) - (95) apply through the 
(n-l)st iterative step, they also apply to the n-th step. Since 
Eqs. (81) - (95) have been shown to apply through the first two orders, 
it follows that these relations must hold for all values of n . 


. t 


V. THE (2n+l) RULE AND SOLUTIONAL ECONOMY. 


The real difficulty of determining the perturbation energies is 
involved in the solution of the partial differential equations which, 
are required. Thus, there is a very real advantage if we can redilce 

/o__\ 

the number of differential equations necessary to obtain e ^2n]f ° r 
e t2n+l]£, ^ rom (2n-l) or 2n to n . This should be possible since 
Wigner showed that if an approximate wavefunction is accurate through 
order X^ , then the corresponding expectation- value for. the Hamiltonian 
should be accurate through order X . Indeed, for non-degenerate 
Rayleigh Schrodinger perturbations, the energy through £ (2n+l) is 
determined by the wavefunction through. ^ ' . We will show in this 

section, how to do the same thing for degenerate perturbations. This 
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solufcional economy is accomplished by means of: (A) The "unwinding 


; (B) ^he repeated use of the magic formula 


and winding" of the 

/ \ 

to "unwind" the > and The sequential solution of a set of 

eigenvalue-eigenvector equations. . 


(A) The "Unwinding and Winding" of the 

3 ^ S J X, 


In carrying out the "unwinding and winding" of the A^ n ^ 0 We 

(p) ’ 

be concerned with two sets of functions: the whicji satisfy the 

(t) 

first s perturbation equations (83); and the l]q w hi c h satisfy 

the first (n-1) perturbation equations (83). Thus, it is obvious that 
if the first n orders of the Hamiltonian are Hermit ean .and - .if u<s or 
if we are dealing with the: u=s is also allowed , then, 

u-1 

1 

P c 


-i 


pf 0 -“lol 11 I[n-l]q UsJ* 1 Un-ljq 


n-u-1 

- I 

v=o 








(127) 


(n) 


The defined by Eq. (92) can be expressed in the form 


£.t 




(128) 


If we apply onto the term with, t=n-l , the Eq. (127) with the u~l 
operation, then we obtain ‘ 


A (n) B Q y 2 y : <(i) (p) ig(n-t-p) . .(t) > 

£l*W . J D p-0 ~ [n-l]q > 


( 129 ) 
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Repeating this operation with u=2 upon the t=n~2 terms in Eq. (129), 


[ •. • < 13 °> 


Thus, after repeating this process (u-1) times. 


A n . - I l «j> p > |S ( ”- t - p >|t t) 

= [s]t t £ 0 p f 0 ~f n_ 


[n-l]q 


where 


u<s for <£r P it or u<s for <i>f p ^ 

~[s]k — 4-[ s J q -» 


l 


- r (131) 


This "unwinding and winding" of has the advantage that instead 

(t) ■ 

of requiring through the (n-l)-st order, it is now only 

necessary to know it through the (n-U-l)-st order. But it has the dis- 
advantage that now we must know the through the u-th order. 

Uslially it is be6t to let u ra (n-l)/2 , or close to it. 


(B) The "Unwinding" of th'e Energy . 

.(n). 

'q - f 


(n) - = » (n) 

~ ln-l]q * [n-l]q 


Since e'V' , = A '71 + e v 'l ; from Eq.. (131), 


itn-ljq " Jq Jq ^(n-ljqln' ’^n-ljq' 


i ' ’ * ! - i > t 

n— u-1 y 


,(p) intn-t-p) (^(t) 


* 


'( 132 > 


where the primes on the summations mean that the term t=p=0 is 

excluded. By taking U=n-1 , we have demonstrated that is 

: sln-ljq 
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ilerniteati' and proved Eq. (100). If in Eq. (132) we replace n '.by 
2 n and choose u=n-l , we obtain 


(2n) 


(0) 


.CnJij.CO) ^ V. n w*(P> 


St2n-l]q " ^[2n-l)ql H l^2n-l]q > + Jj ^[2n-l]q 


!g C 2 »- t -p>| J CO 


[2n-l]q 
(133) 


> . 


While if we replace n by 2n+l- and choose u=n , we obtain 


(2n+l) „.(0) 


<a"" | H ( 2n +l) i,(0) r, V» <d) (p) in(2n+l-t-p) | , (t) 

~{2n]q < ~ [2n]q ■ - ^njq* + J Q 2n] q ‘ U ^[2n)q > ' 

j . • • . ' • ;• V" ' ' (134) 

• - - . 

These two equations demonstrate the (2n+l) rule for the energies? 

*»- 

The ^[2n-l]q anc ^ ^[2n]q ' ^th . t<n in Eqs. (133) and (134) can be 
expressed in terms of the . with m'<n'<n .. 1 


(C) The Unwinding of the * '* ‘ 

v To demonstrate the (2n+l) rule for the wavef unctions , we note 
that in order to calculate according to the magic formula (91), 

we require the matrices 


(n+j) 


.( 0 ) 


a I -x. CO) >,* v * 11 J/ Y 1 oe\ 

~ [n+j-m+p]k ~[n+j-m+p]k ~ [n+j-m+p-l}q ~ [n+j-m+p-l]q * ' 


(n+j) 


which, according to Eq. (92), involves up to (n+j-1) order wavef unctions. 
However, we may apply Eq. (131) to obtain - 


A (n+ ^ 

~{n+j-m+p-l]q 


n+j-l-u u , . 

Jo Jo ^In+j-nrtT-'lJq 


|H (n+J .: t - s) | }l ( ^... lh > , U36) 
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where u< (n+j -m+p-1) . Thus, we can avoid wavef unctions of higher than 
n-th order if we choose u so that 


(j-2) $ u ^ n 


(137) 


For example, taking ipj-1, 


,(n+j) „ V ~y <6 K8) ,guvt-J-c-s; . u, > 

~ [n+j -m+p-1 ] q t f 0 s ^ 0 ‘ ~[n+j-mfp-l]q 1 . ,'2[n+j-l]q 


5 3 


_ n 


» 


, 5 (n+j-t-s)| .(t) 


(138) 


Of particular interest is the case j e n , which establishes the 

(2n+l) rule for , and the case j=n+l , which establishes the 

~L2n-lJq • ' 

(2n+l) rule for "* ' J' ' !l c 

•- ■ i .-■■■■ a _.\;V 

. ■ 

VI. COMPUTATIONAL PROCEDURES: THE BASIC CYCLE AND THE (2h+l) CYCLE ^ 


• .V 


There are many ways in which our formalism can be used to calculate 
the energy and wavefunctions of degenerate states. The optimum pro- 
cedure depends upon the detailed branching of the energy of these states; 
whether the wavefunctions as well as the energy are required; whether 
the energies of all of the degenerate states are being calculated; and 
to what order should the energies and/or wavefunctions be determined. 

In arranging the computations, consideration should be given to two types 
of procedures: the basic cycle and the (2n+l) cycles. In the basic , 

the determined from a knowledge of the <j> j ^ where 

0£m$n-l . In the (2n+l) , the energy through e [2n+lJ2 t * eterm * n ® < * 

(s) - ' : * 

from a knowledge of the $j2n]q ^ or s £ n • If. the determination of 

,(s) 

requires the solution of a partial differential equation and is 
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therefore a computational bottleneck, then the (2n+l) cycle is 

(s) 

beneficial. If, on the other hand, the is being approximated 


by a. linear combination of basis functions [see Section VII] and R 


( 0 ) 


becomes a finite' matrix, then the advantages of the (2n+l) are not 

, (2n+3 ) 

obvious. And, finally, if the energy through the e [2n+l]? re 1 u: *- re ^! 

:.(h+i).; 


then one might choose to use the basic cycle through 
(m) 


and the $ [n]£ 


[n-fl]& 

and complete the calculations by a ‘modification of the 


(2n+l) procedure. ” ' 

In any case, the <{>'s are related to the wavefunctions by ‘Eq. (82) 




(m) 

[n] l 


,(m) 


.( 0 ) 


• [n]£ 


m-1 


t “ 0 - [n— nri-t)q I [n-nrt-tjq 1 v [n) S, 


( 82 ) 


(m) . " 

where the intermediately normalized. For any other normaliza- 

tion, Eq. (86) should be used in place of Eq. (82). 


The Basic Cycle • 

The basic cycle consists of the following chain: 


.t> ef n ^ , , 4>r°i . for i«jm<n-l--&» ^ 

. ~[n-l]q [n]Jl I[n]x, I[n]l " ' - I.[n]x. . ^[njq. 


The starting sequence is the regular basic . cycle with 


.( 1 ) 


( 1 ) 


,( 0 ) ■' .(1)' *-. -.(2)' ■ (2) . .( 0 ) . ,(1) 
?[0)q^ f.[i]f, » $[l)Jl^;f[l)q r^|[l]q“^ *[2]l '» ${2]*"* $[2]£ 
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The steps in the chain are: 

1) The e[^i] q is 8 iven by Eq. (88), 


(n) 


,( 0 ) 


n-1 


5 (t),(n-t) , „(n) . (0) 


. = [n-ljq iln-llq 1 ^ ~(n-l]q ~[n~l]q f! 


( 88 ) 


2) The £ [nj 9 is given by Ec I* (87), 


(n) (n) | - ‘ 

§[n-l]q £ [n]£ i < . . ° *: 


(87) 


3) The expansion coefficients are 8^ ven as the eigen- 

(n) ' V. ^ 

vectors of e j n _i] q by E< 1« (89), , 


A 


(4“^, - 1] <4>r°\, I <J> r°T 0 > « 0 

l ~[n-l]q [n]£ X [n-ljq 1 X [n)£ 


(89) . 


.( 0 ) 


«,i v. T 


The are then given by Eq. (90), 


_ * , - I <*• 


A ( 0 ) ,(o) <A ( 0 ) 1 ,( 0 ) ' 

~[n)£ ~[n-l]q $[n-l]q‘$[n]Jt * 


(90) 


.00 


4) "The ^ or l&m^n-l are given by Eqs. (91)*- (93) , 


I ,(m) 


where 


LOO . y “y 1 ,(p) £(-n-hn-p); (n) 'L.(O) I, 

L~ [n-l]q 4 ^4^ ^fn-n+pj^ k .~In-m+p]y^[n-l)q'^ 


( 0 ) > 
[n]£ * 

(91) 


A<"> . - f ' • 


~ [ n-m+p ] k ; fn-m+-p]k 1 


(92) 


\V 1 
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and 


S (-n+m-p) _ , (n-m+p) ' (n-m+p) ,-l 

k " 1 q “ E [n-mfp]k J * * 


(93) 


(n) 


5) The is the solution of Eq. (83) 


n-1 


fi (°) A ( n ) + T + fH (n) 

H ~[n)£ + -. t = 1 ^ “ e [n)il ] ^n]£, 


(83) 


together with the orthogonality requirement of Eq. (85), 




(85) 


However, if ^ [ n ] * s being approximated by a linear combination of 
fixed basis functions, then is a matrix [see Eq. (.158) J and it 
is more convenient to determine (ji! n L- "b£ Eq. (94). 

~inji. 


■ The (2n+l) Cycle 

The (2n+l) cycle consists of the following chain: 




[2n-l]i for ^ 


„<2n+l),, ,(m) , ■ , - > (2n) .(0) , • (2c) 

*** ~[2nla ^*(2 n U for . ;L S m ^^ E I2»U > 



[2nJ£ ’ I[2n]A^ ~{2n-l]q 



(2n+l) .(0) ^.(n+1) 

e [2n+l}£ ’ 2[2n+l]i T ?[n+l]q 


The starting sequence is the same as the regular (2n+l) cycle with the 
following exceptions: In the n=0 cycle, only the (8)>-th through 


1 
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(10)-th steps are performed; in the n=l cycle, the (2)-nd and 
steps are omitted; in the n=2 cycle, the (2)-nd step is emitted. 

Figure 1 shows the energies and t^’s which kre calculated in each 
cycle. 

The steps in the (2n+l) cycle are: 

1) The determined in exactly the same manner as in step (5) 

of the basic cycle. 

/ \ 

2) The ^[ n +j]£ 5 or are determined by successive applications 

(starting with j=l ) of the following relation (taking m=n ) which is 
obtained from Eq. (91) with the help of Eqs. (135) and-(138). 


?[n+j n 


4> 


(m) 




^[n+j-l]£ , 

+ T m y 1 6 (p) cf-n-j+m-p). (n+j) 

£ p £ 0 ~[n+j-m+p]k G k #[h+j-m+p]k 


< ~ [n+j -1 ] q ^ ~ (n+j ] £ > » 


,(°) 

tn+; 

(91) 


where 


p (-n-j+m-p) r _ (n+j-m+p) _ _ (n+j -m+p) ,-1 
°-k ' q e [n+j-mfp]k J 


(93) 


- < < f ) CO) 1^(0) ' ‘ » 

* r “ • " Sln+j-nrfrjk'rCn+j-m+p-Dq^-fn+j-m+p-l]!! 


(n+j) 
~[n+j-m+p]k 


(135) 


and ' 

.(n+j) r ^r 1 <.(8) i = (n+j-t-s) | . (t) n- 3 Q\ 

~{n+j-m+p- 13 q t “ Q b “ q $[n+j-nri-p-l]q* ~(n+j-l3q * * 


3) ihe tSiu for l^m^I are determined from the same equations, 
(91), (93), (135), and (138) as occur in step (2) where now. j=n-l . 

4) The e [2n^ijq is determined by Eq. (133), 
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P (2n) 

~[2n-l]q 


<( *.(0) Ih^U^ > 
^I2n-l]q' h *$f2n-l]q 


n n-1, 


+ F T ; <t{S’. 1 ,JB < * B - t -» > ii t (t> 


t=0 p=0 


[2n-l]q 1 


[2n-l]q 


(133) 


Here the prime on the summations indicates that the term . t^Op is 
omitted. 

5) The e [2n]i -*- s ‘ ^ eterm ^ nec ^ by t ^ e same equation (87) as occurs in 
step (2) of the basic cycle, except that now n is replaced by 2n . 

6) The is determined by the same equation (89) as occurs in 

step (3) of the basic cycle, except that now n is replaced by 2n . 

7) The ^ or are determined from the same equations (91), 

(93), (135) and (138) as occur in step (2), where now j=n . 

8) The e j2n]q^ is determined by E 9- (134), 


(2n+l) 




( 0 ) 


(2n+l) i . (0) 


•-.t i.- 


tt2n)q 


^[2n)T 


n 


+ r i^<<i> r ( p ) 1 w|H ( 2 n+ 1 ~ t ‘- p) u ^ ) 1 > . 

t= 0 P=0 ~ f2n]< £ ~t2n]q * 


(134) 


Here the prime on the summations indicates that the term t^O^p is 
omitted. 

9) The e [2n+l)& determined by the same equation (87) as occurs in 

step (2) of the basic cycle, except that now n is replaced. by (2n+l) . 

10) Finally, the determined. by the same equation (89) as 

occurs in step (3) of the basic cyclfe, except that now n is replaced 
by (2n+l) . , ' lz r 
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VII. RESOLVED DEGENERACES AND SPECIAL CASES. 


As the calculation outlined in Section VI proceeds to higher order, 
•the original block of degenerate states is split into smaller blocks 
and eventually becomes non-degenerate (or reaches the degree of degeneracy 
which remains in infinite order) . All of the equations and procedures 
in Section VI remain valid when this occurs (and, in fact, are correct 
even when the degeneracy is resolved in zdroth order, i.e. in the non- 

i # - 

degenerate case). The character of the equations changes somewhat, 

. •' 4 ; . . v _ 

however, and the procedures becomei simpler. .* 1 


The major new feature which emerges when the degeneracy is -resolved 


is that the wavefunctions 1 p 


(n) 


begin to be completely determined. When 


q 1 
the degeneracy of the state q of interest is resolved in r-th order , 

the original block of unperturbed states is resolved into (r+1) sub- 

spaces: $[2)k» ••• $[°]k ’ The subsetB £{2-l)k’ $[2-2Jk’ 

are empty. Then the magic formula (91) gives 



(n) 

tr+n)q 


An) 

^[r+n+s]q 


n^O 


s>l 



‘ (139) 


60 that Eq. (82) gives 


; . ♦&!, . < U0 > 
In order to discuss the simplifications which result in the pro- 

i. ' 

cedures of Section VI, let us first -consider the special cases when the 
degeneracy is resolved in first or second order and then discuss the 
general case. 
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(A) Degeneracy Resolved in First Order . 
In this case Eqs. (140) and (91) give 




(n) 


A (n) 

*W]q 


*< n ) + y *(0) p(-D A (h4-i) 
*[n]q + l ?[l)k C k tU)k ' 


(161) 


where, according to Eqs'. (140) and (94), 


(n) R (0) j j(t)^(n-t) 

~tn]q t“l • q 


4> 


(142) 


and, according to Eqs. (140) and (92), 


im” * ltS.> + T <f g} JS (t> • . WJ) 


s[l]k 


:[i]k' 


[n] < t=2 2 t« k . 


The computational procedure involves solving Eq. (142) for 

which allows to be calculated. Then the (2n)^th and (2n+l)-st 

Q 

V 

order energies are calculated by Eqs. (9) and, (10). 


(B) Degeneracy Resolved in Second Order . 
In this case Eqs. (140) and (91) give 




a) „ ad „ ad 

q ~(3]q ~t 2 ]q 


k 


(-2) .(3); 
~[2]k 


and 


A.a> 


a) 

[i]q 




(0) £(-l) (2) *? ,t.u; i, 

X[ 1 ]A tniHl ^[l]q‘^ 


.(0) u(0) 


U2]l 


(144) 


(145) 



1 


Then the higher, order wavefunctions are given by 


|.(n) ,(n) = , (n) r,(0) e(-2),(n+2) 

q ~[n+2]q ^[n+l]q + £$[2]A ~[2]k ’ 


( 146 ) 


. (n) 
•[n+l]q 


,(n) y,(0) c(-l) (n+1) y,(l)=£(-2) (n+1) 

?[n]q + itl]k + [2]k K • £[2]k 


(147) 


where the coefficients in Eqs. (144), (145), (146) and (147) are 




‘ft? - ' - ’ a48) 


$ 5 ? • 


The energies are determined by Eqs. (9) and (10) 


(C) Degeneracy Resolved in r°th Order . 

In the basic cycle, if the degeneracy is resolved in r-th order we 
immediately obtain hy Eq. (140). Then in each succeeding cycle, 

we obtain one additional perturbed wavefuncticn. The only modification 
of the procedure is the omission of step (3). Of course, steps (1) and 
(2) are much simpler because e [^i]q I s a nuIfl ber, not a matrix. In 
step .(4), it is necessary to consider only n-r$m£n-l and the upper 
limit in the sum over p in Eq. (91). becomes m-(n-r) . ' ' - • ' ’ 

In the (2n+l) cycle, if the degeneracy is resolved in r-th order, 
again we immediately learn • In the next cycle we obtain 
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i 

) 

. . In . ' 

and . ^ * ; and/feach subsequent cycle, two additional until 

I q 9 

n*=r . Thereafter, each additional cycle yields one additional 
The modif ications in the procedure are that steps (6) and (10) are 
.omitted, while in steps (3) and (7), it is necessary to consider only 
those perturbed functions which are not complete according to Eq. (140). 
Thus after . is obtained, steps (3) and (7) are omitted .entirely. 

v 

When the perturbed wavef unctions are known completely, Eqs. (133) and 

(134) reduce to Eqs. (9) and (10). 

* / 

The computational procedure fhr the case .when the degeneracy is 

i 

resolved in fourth order is illustrated in ; Figure 1. w 

I • . *■>' 

> . - *T . 

, ' . - ' * - Y ' 

• * . • ,'> . . ' - ^ 

• .... . •• .. X - 



Figure X. Computational Procedure . This figure illustrates the quantitites obtained in each 
computational cycle. The basic cycle is indicated by the dotted lines, while the (2n 4- 1) 
cycle is indicated by solid lines. See Section VI of the text. The. illustration is for the 
case when the degeneracy is resolved in fourth order. 
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1 


VIII. VARIATIONAL METHODS FOR APPROXIMATING THE <J> 


(n) . 

[n]q * 


As mentioned in Section V, the most intractable aspect of degenerate 

perturbation problems is the solution of the partial differential 

equation (83) with the Eq. (85) constraint (or equivalently, the solution 

/ \ 

of Eq. (94)) for the function • The problems involved are similar 

to' those of determining the n-th order wavefunction in a non-degenerate 
example. In this section, we consider variational methods of approximation 
which can be used when Eq. (83) cannot be solved exactly. In this section, 
tildes are used to indicate approximate functions or energies. -v 
The basis of the variational method is the usual variational 

*• . "Y ’ 

principle: Define the functional J as . . , , . 


J - <$ q |H(A) - E q (A)|ij} q > . , 


(149); • 


«. „ .. M-.t . - ■ . ‘ 

Holding E^ fixed, vary so as to make 6J = 0 . Repeat this process 

for different values of E q until simultaneously 


J *= 0 and <$J *=> 0 .j 


(150) 


The application of variational methods to degenerate .problems has been 

12 

considered previously by Epstein. * . 


(A) Linear basis set.' 


A generally applicable method is to express ty q as a linear 

combination of the and a fixed set of basis functions f, . It 

~ [OJq k 

is convenient to express the set of f^ in vector form as f . Without 
loss of generality we can require* that the f be ortho-normal, orthogonal 
to the > and diagonalize , 


54 


<f|f> 




(151) 


and 


<f k |H< 0 >| f .> - £<\ 


(152) 


The effect of the linear basis approximation is to reduce the infinite 
dimensional Hilbert space to the sum of the two finite subspaces, f and 
^[0]q * P ro 3 ector onto this reduced subspace is 


P “ l£ >< £l + 


(153) 


In this notation. 






£<£\l> + <j> r ^«J> 


~ ~ q 


<°> <*<« |J 

~tO]q i[0]q ,v q 


li* ~ 


(154) 


Here the expansion coefficients <f|^> and << ^[o]q^q > are var l at i° na l 
parameters. Substituting Eq. (154) into Eq. (149) and satisfying the 
variational condition Eq. (150), we obtain the familiar eigenvalue-eigen- 
vector equation 


[H(X) - 


E 

q q 


(155) 


where 


H(X) 


PH(A)P and 


H^ n) = PH^P 


(156) 


From Eqs. (151) and (152), the zeroth order equation becomes 


[H (0) - e£ 0) U k = 0 and [5 (0 > - e< 0) ]$[°] q = 0.(157) 
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The resolvent corresponding to H v is then 


l ' 4 0) - • 


(158) 


Thus, if PF is an arbitrary function in the reduced space then 


-(H <0) - e< 0) >£ <0) PF - 


(159) 


The resolvent R v then plays the same role in the reduced space as the 
( 0 ) - 

resolvent R plays in the full Hilbert space. 

Thus, we irfay apply the formalism of Sections II- VII to solve JEq. (155) 
exactly to any desired order by using the tilde matrices in place. of the 
corresponding operators and by using PF in place of all functions F . 

i 

Thus, the solution of the degenerate perturbation equations is reduced to 
the solution of a set of matrix equations. 

: • "-.t v. !“ * 

( 2 ) 

(B) Variational principle for e [i]q * 

The general application of variational methods to degenerate problems 
leads to considerable ambiguity. Here we restrict our consideration to 
the determination of the energy through second order. 

Let us consider a trial wavef unction of the form’ 


? -• + X* (1) , 

q q q 


(160) 


where, as suggested by Eq. (28), 


,( 0 ) 

hllq^Ulq 1 ^ > 


( 161 ) 


and 


1 


,7,(1) _ I<1) <A (0) |,7<0). , .(0) (0) ,7(1) 

^q " ~fl]q ~U3q' V q ^[0]q^[0jq^q '■ . (162) 

The notation implies that we know the sets ^ [ 0 ] q an< ^ ^[l]q exac tly» 
and the quantities with tildes are to be obtained from the variational 
principle. Substituting Eq. (162) into J and expanding in powers of 
X gives as the leading term 

J«> - <$ (0) lC, -? <2) )<4?> |i <0) >; 

r q ~(lJq ~[i]q q ~ Ci3q 1 q 


where 






[ 1 3 q T ~[l]q 


~ [ilq 


(163) 


Note that when is exact, becomes exactly. Also, 

~[l]q ~ 1 1] q •• -t . - ~[l]q 


the J 


does not involve the coefficients 


<♦£> li (1) > 

"[0]q^q 


at all. 


The variational equation for one of the members of the set » 

say , is obtained by putting <$J^=0 for variations 

4 *[l]qa ' ThlS glVeS 

' * <2) ■ ' £*sl is (o >us > + >>*fs3 q nr* 


+ complex conjugate = 0 


(164) 


The variational equation for the <t f > [l]qJ^q^ > t ^ ien obtained by 


variations ^ << ^[l]ql^q^ :> » which give 


■l«g> |* (0) > 

a [l]q q Hl]q'-% 


( 165 ) 
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Thus, as discussed by Epstein, the equations for the a ^d t ^ ie 

< *{Sq | C > are coupled together and presumably must be solved by an 
iterative procedure. 


(c) Quasi-variational procedure . 

It is easy to write down a set of functionals whose stationary points 


/ \ 

are the exact solutions <b r ■, 

~lnjq 


(n) 

[n]qa 


If 3 mem her of the set, then 


j(2n) 

■ qa 


+ X {< *Wq«l s<t> -l*{SS > + c -°- 1 a66) 


is such a functional, where we assume that the lower order functions are 

/' 2 n \ 

known exactly. The most straightforward procedure is to use _ <$J V =0 to 


1 60 


qa 


determine ^ [n ] qot * an< * t ^ ien to us3 this set in the computational procedure 
outlined in Sections VI and VII, treating it as though it were exact . It 
should be emphasized, however, that this.i£ not a truly variational 
procedure in the sense of making the total energy stationary. 


IX. THE OPERATOR TREATMENT. 


29 


In the present section, we wish to develop the operator approach 

f 

independently of the direct treatment and then show how the two are related 
A. Definitions and basic strategy . 

Before proceeding with the derivations, let us introduce some of the 
new Hermitean operators which we: require: First of all, there is a 

family of resolvents .■ ' * ‘ 


* ( - s) 1 %-’ 


(167) 
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Then, letting V » H - , 


l [(V - E n + e n (0) )R (0) ] t (V - E + e (0) ) = l X°Q^ n) . 


t=0 


q q 


q q 


n=l 


Here the n-th order of is given by the recursion relation 


Q<"> - iW .+ "I . 

U t=l ' 0 


(169) 


And, finally, there is the set of operators 

% : K P + X “V* - - V ■ X A P n> *, 

• ' . * i w 

where ' 


K 1 = and K p 


Q O l)+ j 2 <-2 • : <m) 


One of the most important properties of the Q q and Q p is that 


% " < v - E q + e q 0) 3tl + R (0) Q 0 3 


(172) 


and 


<L , + (Q , - K )R ( " P) (Q - K ) . 
Pi P~1 P P P •» 


(173) 


M): 


Using Eq. (173), it is easy to show by induction that the Q p ' are 
given by the recursion relations, 


Q P 


(n) 


VI 




n<p+2 


& 


‘ <£° ' - + T Q'!^ ( ' f) q' n+P ' s) . 

* * R=n+1 ^ - ■ 


+ qi-3 1>R< ' n+2>< 4-3 1> 


n-1 

I 

s=p+l 


n=p+2 

n>p+2 


(174) 


(168) 


(170) 
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In order to explain the significance of these operators, let us 

• r • • • 

briefly sketch the strategy of the operator development which is given In 


subsection B. Our attention is focused on the wavef unction , which we 

relate by means of infinite order operators to the unperturbed functions 

. The most general form for lb is 
~[0)q & . q- * 


« - R<w V*[S], < i[S),IV ’■ 


(175) 


where Qq is an operator to be determined so that satisfies the 

Schrodinger equation. Here might be any of the functions ^[o]q * 

Since annihilates all of the functions $[o]q » we com P^- ete the 


definition of Qq by 

. x 

^[Ojql’oUlSlq^tOjqlV ‘ 0 " ^ 

Thus Qq becomes an effective Hamiltonian -operator in the space ^[o]q * 
This is not a true secular equation, however, since Qq is not linear in 
E^ . However if we consider the leading term in Eq. (176), we obtain the 
standard eigenvalue equation 


<A0) | n (D|,(0) .(0) 1.(0). _ 

$[0]q' Q 0 ~ [0]q >< ~ [0]q* ^q 0 * : (177) 


H (1) - £ (1 > 


(178) 


This is identical to Eq. (37) of the direct treatment and defines the 


resolution of 4 rn i 

. - ~[03q 


• -(b) (o) 

onto the subspaces and ^[l]k ’ ^en we 


resolve ^ ' in the more specific manner. 


so 


- a + R^Q^a + R^Q^t^tfJjqlv 


(179) 


where is an operator to be determined so that satisfies 

Eq. (176). Here \j> might be any of the functions this 

way ivc proceed progressively to resolve the degeneracy in exactly' the 

same manner as in the direct treatment so that 'J> r 1n and the 

Injit — l n j it-. 

have the same significance in the two formulations. If the degeneracy is 


resolved in the r-th order, the final expression for \|> is 

• ip * n a + r ( " p) q )\p (0) , • 

* P =o . p « 


(180) 


(0) 


where the is an effective Hamiltonian in the subspace ^[p]q * 


’ 0 - 


(181) 


A key step in the development is to show t'tjat . is given by Eq. (170) 

or (173) and that, as a result, the expression for can be 

rewritten in the more useful form of Eq. (231) below. 

From this analysis it becomes evident that 


bf % * 
(182) 


so that ef n |. and* 4>r 0 io are determined in exactly ‘the same manner in 
[n]5l ~[nj& 

the operator and in the direct approaches . The operator equivalent of . 

, (m) , 

9> r •>„ is obtained in subsection C. _ 

~ [njS. 

We turn now to the detailed derivation of the operator treatment. 

By Infinite Order -Perturbation . Theory ♦ • , . 

. . ‘ * • .1 

We start with the unperturbed Schrodi'nger equation 


(H ' E q ^[0], 


(183) 
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At each stage of the development, we seek the energy and the wavef unction 
I p corresponding to the perturbed Schrodinger equation 

q 


(H - 


am 


The most general functional form for the wavefunction is - A -.. ; 


fp (0) T 4- n A (°> ^(°> |, (1 ' 

tR J 0 + ’ 


(185) 


where Jq is an operator which we wish. to -determine. Since ,R V ' is 
orthogonal to the space of the ^[o]q * without loss “Of generality^ ve 
can require that *:”• . ■ 

^[0], I v. ■. 011 '■■■■ : <186) 

(0) - - ■ Co) (0) 

Since «i> r _ , . \can be resolved into the spaces of d> r \ apd : £, <j> r 

~ [Ojq r ~[n]q s=l ~[s]k :/ - . 

Eq. (186) is equivalent to the set of equations 3 ) 


<4>r°l l J nl<f>rnl U > " 0 

~ln]q ' 0 ~[0]q ^ £ 0 ] q * ^q 


(187)’ 




l<s<n . . 


(188) 


Applying the operator onto ijr as given by Eq. (185) 

fi(0) ^, * S<0,E<0 %i[S]q<*[?fql«:q> " - <189 > 

* ' \ _ * . 

But rewriting the perturbed. Schrodinger equation. 


O -[V - E + 
q q q . r q 


- t v : E q + e q C0) HR (0) i- ; Uf[S] q <j[oj q IV • 


.(190) 


1 


Then, comparing Eqs. (189) and (190) ,' 


J 0$h < *[0]q!’C )> - 17 - E o + E a <0)l!1 + E< ° )j 0 I ?(S]a^fSi a IV ' tl91 > 


q q 


o J S[03q"S[o] q 'V 


Applying the operator on both sides of Eq.(191), we get a Lippmann- 


„ . 30 

Schwinger type of equation , 


a = b(c + a) “ b[c + b(c + a)] = b[c + b{c + b(c + a)}] * 

- . CO ; 

' s ... I (b) t c .. . •' (192) 

• . ’ t=l 


Thus, we obtain 


; . ' R< ° )j 0 ![ 0 ] q < ?[ 0 )IV (193> 

where Q n is the infinite order Hermitean operator given by Eq. (168). 

U and • 

Substituting Eq. (193) into Eq. (191)/ making, use of Eq. (172) , 


JL<b^ <cb rn , U > Q„d>^ lip > .« 

0i[0]q ~ [0]q ’ *q V 0M0]q i[0]q |V q 


0~[03q ~ [0]q * *q 


The wavefunction (185) can then be written in the form 


(194) 


V [R<< \ + ^[OJq^lOJqlV . a95) 


i Also, the condition on Jq becomes 


< ![03ql Q 0^[0]q >< ^[0]ql V = ° 


( 196 ) 


The energy is one of the roots of the "generalized secular 

equation” ", . - 


< ~ [ 0] q ^ Q 0 ^ [ 0 j q > I ' “ ° \ 
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Eq. (197) is not a true secular equation since Qq is not linear in . 
Except for notation, Eqs. (195) -(197) have been derived and rederived many 
times in connection with "infinite order perturbation theory". If 
were non-degenerate, then ** ^ [ 0 ] q * ^(q) = * ,an ^ 

«* <rp W \Q^ + • For degenerate states, it is not 

99 9 9 resolve the- 

necessary, but it is convenient, to successively/set of degenerate states. 

(1) First Resolution . 

Iiv the" limit that X becomes small, Qq becomes Qq^ “ and 

ip becobes . In this limit, Eq. (196) becomes 


«j><0> ><*<« U<°>> - 0 . 

~ 1 0 3 q ”0 . ~ [Ojq ~ [ 0 ] q *q U ’ 


(198) 


Of course we recognize that 


flO]q 


..(0) J.(l) (l) u (0) . nQQ , 

< ~[0]q Q 0 + .?o,.J-$r 0la > (199) 


<T- • >. ~[0]q 


i and Eq. (199) is exactly the same eigenvalue-eigenvector equation as the 
i first order Eq. (99). Thus, in exactly the same manner as in the* direct 
treatment, we determine the anc * t ^ ie corresponding • The 

! adjoint of the eigenvalue-eigenvector equation Eq. (102) is then 


a (D 

sll]*- 


- • < iujf |, o 1) IJ[S], > * C - e o 1>)< triu If roi/ <200) 


UU ~[0]q~ 


However, when s > t 


(0) - a<°> <*(0) , .1(0) ><A (0) I (0) VO) ,,(0) . 

-(t3q < -[t]q ~[t+l]q ^tt+l]q' ' * ‘ ^[s-l]q ^[s-l]q ^[s]^ * (201) 


X v 

~[s]£ 


It follows that if we "back-load" Eq. (200) by using, Eq. (201) when 
s > 1 


> 
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.*& (202) 

Then, if we "front-load" Eq. (202) by using Eq. (201) and Eq. (84), 

‘ <e [W- E “ )S .,i i t ,A,ri' s>i, t>i . (203) 


Eq. (203) can then be stated compactly in the form of two relations 


R( " S)< 2o 1) ![t]q “ ° » t l 1 


(204) 


R (-s) Q (l) R (-t) 


-R (-1) <S _<S - , s>l, t>l /■ - (205) 

S j !• ji X «i . 


Now let us return to the problem of completing tfie first resolution. 


If we define R v X 'J 1 so that 




(206) 


Then, resolving the ^[o]q onto the s P ace the > Eq. (1$4) 


becomes 


MiWSh'V ■ + • ; , < 2 ° 7 > 


[ijq ariUq'V 


It follows from Eqs. (188) and (207) that 


[* 4 ' w <q 0 - q‘ x) ) + r^q^ux + <208) 


(-i) T ,,(0) ^(0) 


Making use of Eqs. (204) and (205) and rearranging, Eq. (208) becomes 

R ' :1 \*uVuhlV - ; eW V- < 1> >t l + R( " 1>J i 1 Sm q < *mqlV (209 > 


1 


• • --30 • •• 

Eq. (203) is another equation of the Lippmann-Schwinger type which has 

as its solution 

• V ■ : 


x< ~X*ujX°) J V ' . " (210) 


Making use of Eqs. (174) and (204), Eq. (207) becomes 

: J oJ[S]< l < +tS],fV " (1 + C^Wulq^UJqlV - (2U) 

Thus, making use of Eqs. (206), (210), and (211) , the vavefunction (185) 

y - T # • * 

V 

becomes ' 


- IR <0) (1 + Q^V" 1 ^ + 1 + ^K^Iljq^UJqlV'' 1 (2l2) 

Also, making use of Eq. (204), the condition (187) on gives the new 

eigenvalue-eigenvector equation for E^ , 


< iu], lQ i | Juj» >< i{u, l V ■ 01 


(213) 


(2) Second. Resolution . . 

K . The subsequent resolutions follow rather closely* to the pattern 
established by the first resolution so that it is not necessary; to give 
as much detail. V * *' 

In the limit that A is small, Eq. (213) becomes 

• ' • _ ' ; ■ *. * 1. 

' ' +x2 ^U], lQ o 2) l*UJ,^*U], l ^° )> +0<x3> ‘ °‘ <2W) 
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According to Eq. (203) the first two terms in Eq. (214) vanish so that in 
the limit that X is small. 


ln (2 >U (0 > ><A (0 > U<°>> 

^U)q lQ 0 'Stljq 


0 . 


(215) 


Here 


q (2) = H (1) R (°) H (D + h ( 2 ^ . Thus, ve recognize that 


e < 2 > = <A<°> !q ( 2 > + e ( 2 ) U (0) > 

zll]q $tl]q ,Q 0 + E q '2[l]q 


(216) 


and Eq. (215) is exactly the same as the second order eigenvalue-eigenvector 
Eq. (99) : . Thus, in exactly the same manner as the direct treatment , we 
determine the e £ 2 ] the corres P OIU ^ n S ^[2]£. * ^- so » ih - - the same 
manner as Eqs. (200) -(205) we can prove that 


•d(-s) (2) , (0) 

R Q 0 ~ [ t)q 


0 , e>2, t>2 


(217) 


and 


••'it 


R (- S ) Q (2) R C- t ) . -R ( - 2) 6 s _ 2 6 tj2 , si2, t>2 . 


(218) 


X-2) 


In order to complete the second resolution, let us define R J 2 
so that 


’ k 


? Juik^UJklV " R< ‘ 2)j 2*mq < *[2kl’ f q > • < 219 -> 


2~[2]q~[2]q l y q 


Then, resolving the space of the onto the spaces of the ^[2]£ * 

Eq. (211) becomes 


V^q^tOlq'V ■ « + ^« )V1 + V (220) 

Thus, using Eqs. (188), (205), and (218), 


i + 


_.(o> 
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. . ; [r ( - 2) (q 1 -k 2 xi + r ( - 2) j 2 ) - R ( ' 2 ) J 2 ]*[ 2 ],<J[ 23 ,IV " 0- <221) 

30 

Eq. (221) is a Lippmann-Schwinger type equation whose solution.,is 




(- 2 ) 


(0) „ A (0) 




(222) 


Then with the help of Eqs. (173),' (204), and (217), Eq. (220) becomes 

” (1 + + V < ' 2> )<>2*ui 4 ^S 4 | V (223) 


Making use of Eqs. (173) and (217), Eq. (212) becomes 


v ' 1 ! 




R (0) (l + + K 2 R ( 2 ) )Q 2 

+ (1 + R <_1) Q 1 )(1 + R ( “ 2) Q 2 ) 


£ 


*{S|q < *Wql*q > ' <22 ^ 




Also, making use of Eqs. (204) and (218), the condition (187) on 
gives the new eigenvalue-eigenvector equation for ,. y - , - *>1 % i '■ 


^mql^'imq^mqlV - 0 • 


.(0) s,A<0) 


<225) 


(3) . The General (n+l)-st Resolution . 

The general pattern for the resolutions seems to be well established 
in the first and second resolutions.' We wish to prove by induction that 
it applies in general. The proof consists in demonstrating that if the 
pattern applies to the first n resolutions, it also applies to the 
(n+l)st . The pattern may be outlined as follows: - 
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First, is the n-th order energy-eigenvector relation 




(226) 


Comparing Eq. (226) with Eq. (99), . it appears that 


e (n) 

f[n-ljq 




(227) 


The real proof of Eq= (227) for a general value of n is made in the subsec- 
tion C. Thus, in exactly the same manner as in the direct treatment, the 
and the corresponding are determined. Also, in the same 

manner as Eqs. (200) -(205) it follows that " . • 


&<e M-k[t] = °* s i. n > c i n ( 


(228) 


E (- S )^n) R <-t) 


/ \ « • - 4 t V ' ?\ 

— r “ 6. , s>n, t>n 

s,n t,n - — 


(229) 


After the n-th resolution, we assert that 


J 0 ![0), < f[0],l*q > ■ (1 + KR ( . (230) 




R (0) { n (1 + K R (-P ^)}Q 
P”1 . P n 

+ n (i + r^“ p) q ) 

_ P =1 p 


C^lSqV 


(231) 




( 232 ) 


This completes the pattern . 
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Now we will .start on the inductive proof. According to Eq. (174), 

(233) 


n-f2 co 

0 - AQ' 11 ' + y A V P + l A p Q< p) . 

^ 0 p=2 P- 2 p-n+3 n 


Thus, Eq. (232) can be expanded in powers of X to obtain 


t=l s=0 


Here, for convenience, we have let = Qq^ • Then, making use of 


Eqs. (201) and (226), 


~ [n]q ^s-2 ~ [n]q 


~ [s]q ~ [s-l]q ~ [s-l]q ^s-2 ~ [s-ljq J 
■ 0 • *- • - ; < 235 > 


A ' 


Thus, the only term in Eq. (234) which has not previously been shown to 
be zero corresponds to t *» n+1 *» s . This means that in the limit of 

. ■ ”-‘t v - ■ . 

X becomes small, Eq. (232) becomes 


<c q iei i) i*s q ><c,x o) > - ° 


h T4 (236) 


From Eq. (236), the sequence of Eqs. (227)-(229) with n replaced by 
(n+1) follows easily. . . 

r V s 

In order to complete our (n+l)st resolution, let us define 

n , so that - : - 

n+1 ; 


. i • <237 > 

Then, resolving the space of the onto the subspaces of the 

~lnjq 1 

(0) .'■ = *■ ' i ' 

$[n+l]£ ’ Eq * .( 23 °) becomes * •• 


.(-n-Dr ;co> ^(o) 
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T A (®) J, x CO) I , > ,. = 

0^[0]q $[0)q IV q 

[ S (1 + K R (_P) ))Q n (l + • (238) 

P = 1 


But according to Eq. (188), if R^ n ^ is applied to Eq. (238), both 
sides of the resulting equation vanish- Making use of Eq. (229), the 
resulting equation is . .... 


+ R< " n " 1)j n + i ) ?[n+l] q ^["l lq l V ' 0 ' 


First, we replace the Q by (Q - Q^ n ^^) + qC 11 *-*-) 

n n n~JL n~l 


(239) 


Then remembering 


that by the Eq. (229) with n replaced by (n-KL) [which we h'avq, just 
shown], R C-n-l) Q (n+l) R (-n-l) = _ R (-n-l) ^ Eq> (2 39) becomes ~ 


[R (_n “ 1> (Q n - K n+1 ) (l + R ( “ n 1) J n+1 ) " R ( ~ n “ 3 ) 


. .. M-.t 


J n+1 ] $ [n+1 ] q < ^ [Si] q ^ V 


(240) 


Actually, in Eq. (240) we have done one more thing. We have replaced 

0 ( n +l) f n tbe f£ rs t term by K . Our justification is that the extra 
'n-1 n+1 J 


terns, K n+ i *" “ K n * not contri ^ ute to Eq.(240) because of 

Eqs. (228) and (229). The Eq. (240) is of the Lippmann-Schwinger type 
and has as its solution 


30 


(-n-1) 


.( 0 ) 


( 0 ) 


R ' " •' J »HliSl] q < *iS-lJ q IV * R< '°' 1) (Q n+l - W+t^D^lSiDqlV 


.(0) ,a(0) 


EC ' n ' 1> V 1 +tSa) q < ^[nii) q IV < 24l > 


Next, let us note that as a result of Eqst (173) and (228) 
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V 1 + R< ~ n ~ 1>( Vi - ‘W 1 * 


(0) 

[n+l]q 


- tu + W W) >Vi : W C ' , \i l !®i], 


" (1 + W[£l], ' 


(242) 


Thus, substituting (the first form of). Eq. (241) into Eq.' (238), we 
obtain Eq. (230) with n replaced by (n+1) . 

Then, comparing’ Eq. (185) with Eq. (231) and using the second form 


j of Eq. (241), 


, ~ [0]q l[0]q ,v q 


n 

n 

(1 + R ( " p) 

p=i 


n 

n 

(1 + R ( " p) 

p^i 


n+1 

n 

(1 + R (_P ^ 

p=i 



V^nJq^lnJq'V 


...t 


n+1 ~[n+l]q ~[n+l]q' r q 

[nllk'V •' <243) 

Thus, substituting Eqw (230) with n replaced by (n+1) together fwith 


1 Eq. (243) into Eq. (185), we obtain Eq. (231) with n replaced by (n+1) . 

} 

* 

; And finally, using Eqs. (228) and (230) with n replaced by (n+1)., 

* * 

* f • 

Eq. (187) with n replaced by (n+1) becomes Eq. (232) with n replaced 
; by (n+1) . This completes our proof! 

We have shown that if Eqs. (226)-(232) are valid for the first n 
, resolutions, they are also valid for the (n+l)-st resolution. And 
since these equations apply to the first and second resolutions, it follows 
that they are valid for all values of n .* ' 


1 
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C. The Operator Expressions for the r™ ^ « and the Operators 

~ [n J a/ 'n-2 

Having completed our formal development of the infinite order pertur- 
bations, let us make a closer examination of the operators which occur in 
that treatment. First of all, we can write Qq^ in the form 


j 5 <a) R C0> B (b>_ R (0) s (c) _ ny ( 


(244) 


r where the summation is taken over all possible terms such that 


a + b + ... + c = n 


i and 


■ 1 ^ a, b , . . * , c ^ n. i 


(245) 


/_ \ 

Also, Q 1 can be expressed in the compact form 
P 




n>p+2 . 


(246) 


Here the summation is taken over all possible terms such that 


a + (b -- p) -+* i — (c— =-_p.)~ 


p + 1 a, b, ... , c n 


(247) 


Another useful expression for is. 


< 


1 


“ H (1) [l + R (0) Q 0 n {1 + R ( ' s) (Q e - K s )}] 

s=r & 


+ l H (t) [l + R (0) Q ] n {1 + R ( ~ s) (Q - K )} 

t=2 8=1 S s 


+ l Q^hl n {l + R (-s) (Q - K )}] 

t*»2 * Z s=t S S 


(248) 


Eq. (248) can be proved by Induction in the following manner: first, 

expand (V - E + e ) in Eq. (172). Then use -Eq.(173) together with 

q q ■ 

the new Eq. (172) to express in the form of Eq. (248). Next, express 

Q by Eq. (173) in which Q is given by Eq. (248)... Finally rearrange 
the terms so that Q is given by Eq. (248). 

/ \ 

Of all of these operators, the special family of Q n _2 is especially 

/ \ 

important. And, of all the expansions of tlie most useful one is 




y h^^g^ i 

t=0 n 


t"? 1 G (t)tjj(n-t) 
t=0 n 


(249) 


Although the H v , the Q q , and’ the resolvents are all Hermitean, the 
operator G ^ is not Hermitean. However, its adjoint is easily 

obtained by reversing the order of the Hermitean operators which form 
G^ fc ^ . The operator expressions for can be obtained by expanding 

Eq. (248). The expressions for G^*^ through n=5 are given in Table 1. 
For all values of n , .. . 


( 250 ) 


For values of n > 2 , we found two equivalent recursion relations for. 
determing the G 1, ' : " " 
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( t) 

Table 1. The Operators through n = 5 . 


.( 0 ) 


n 


= 1 


for all values of n 


.(1} 


(OKU) 


R' Q 


0 


= R^Q^ + r (_ 1 ) q^ 2) 

S< 2 > - R^U + Q^R^lQf- 


,( 1 ) _• 
J 4 

,< 2 ) 


R (0) (1) + r (-1) q C 2) + R (-2) O) .. 
x 0 0 U 


- R 


(0) [( 


l + Qo R 


UK(-1) q ( 2) + Q CD R (-2) Q (3)j 


,<3) 


+ R^U + Qq 2) R (_2) ]Q^ 3) 

= R (0) [l + Q^V^m + Q^ 2) R ( : 2) t ]Q.[ 3) 


r (l) 

G 5 


.(2) 


- R (0) Q^ 1) + R ( - X) Qf + R ( - 2) qf > + 

« R <°>[q< 2 > + Q< 1 >(R ( - 1 >Q< 2) + R ( “ 2 > Q f + R^qfb] 

+ W(i + qJ 2 T r" ( ^ 3 ^ 2 T r^ 4 T ; ' 


3 i 


+ e ( - 2 >U + q' 3) s ( - 3) Jq' 4) 


J ft lf “ 


.(3) 


R 


( 0 ) 


(1 + q£ L) R ( ~ 1) )[(1 + Qg 2 ) R ( ~ 2 ) )Q^ 3) + Qq 2) R (-3) q| A: 13 


+ Q^ 1 ) R C 2) ll +'q< 3 ) R ( ^ 3 ) ]q< 4) 


;+ .R ( - 1 )u + Qf>R ( - 2) ]U + Q< 3 >R ( ^]Q^> 


.(4) 


R (0) [l + Qq^R^ ^ni + Qq 2) R (_2) H1 + q^R^-lQ^ 


1 


(1) The first set of relations is obtained directly by using Eq. 

(2) (3) 

(248) to successively expand ,Q^ »••• • In using Eq. (248) to 

.express Q ^9 * °nly the previously expanded are needed for the 

right-hand side of Eq. (248). Thus we obtain 


R (0) Q a) + l , 

S=1 


,(-s)„(s+l) 


. |G (t ; 1) E ( ‘ n+J) 
n-2 


- G (t) . G (t > + 


+ K " * + - 


p=n-s 


2 < t < n-2 , 


G (n 1) ■- R W H (1 +' ( Q (p ^R ( " p) )Q (n “ 1) . ' (25: 

n P = 1 P~** n-3 

.(2) The second set of relations is the operator equivalent of the 

... . .. , r .. '. : •' M ! t -: 

magic formula 


,(n-l) 


r.®> «■ T 

L p=i 


p+1 . J H n-3 


+ l G (t " p) R 
n-1 **, n-p 


(t-p) (-D+l+p) (n-l) o < t < n-1. 

n-p n-3 


Eqs. (254) and (255) can be proved by starting with Eqs. (251)-(253), 


making use of the identity 
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n-2 

n 

p°n-s 


(i + c< p _y- p) ) i + <>ir s s y~ n+s) v 


+ °f "T 3 <1 + Q W R <-'’> )q (n-2-u> E (- n +2+u) ( 

u=0 V«=n-6 


(256) 


interchanging the order of the summations over s and u , and then 
making use of Eqs. (251) to (253) to identify the terms which go together 
to form G’s on the right-hand side of the resulting equation. 

The point of the greatest interest is that i 


(m) ■ h r (m).( 0) 

?[n)£ . G n+l^n]E * 




(257) 


»•* 

Eqs. (250), (254), and (255) provide the operator equivalent., of the 
magic formula. On account of Eq. (256), Eq, (249) gives 


n 


■V .. i H <c> 4 n -F> - q (n h<°>,, 


(258) 


.Thus, 


i (n) _ ■ .vu; i , (.u; • 

tw $[s]* lQ n-2 ! ^[n-l]q * 


.( 0 ) 


(n). , (0) 




(259) 


and our whole direct formulation can be transcribed into the formal 

• .... 4 ■ > (n^) ■ * ' * * - 

operator language. The "unwinding" of the i 0 as given by Eq', (90) 

z I s i*' 

corresponds to tfce unlikely-looking identity - 


u n-l-u 


0 c?> c 

V *2 * L n ** 8+1 

p*0 :t»0 


r i g^w, 

_ n 8+1 U 


(260) 


a(h) 


where u < 6+1 < n , If s > u , then the expressions for are 

n*”4 

more complicated than the ones obtained with a = u , However, the . 
additional terms cancel each other because" of Eq. (205). 



1 
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In the order-by-order development of the energies and the wave- 

‘ . functions, it is the functions G^<j>^io = which are determined 

n+lL[njx. n]£ 


The 


and these functions, in turn, determine the . 

n-1 Lin]/, 

determination of the ^n+l^[n]£ anc * t * ie order-by-order development of 
the operator formalism proceeds in exactly the same manner as for the 
direct formalism given in Sections III-VII. If r is the order in which 
the degeneracy is resolved, then . Q n = Q r if n > r . Thus, the 
relationships between our present work and the. previous formal; operator 

* * ' : * . " i 

treatments becomes apparent. The. one approach can be completely transcribed 
into the other. ' . . 
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24 Throughout this paper, w a use the convention that whenever the upper 
hound of a summation is less than the lower bound, the sum is equal 


,<0> 


to zero. 

25. If the Hilbert space is finite sc that R''*' is a known matrix, 

r^ D ''f can be obtained simply by matrix multiplication [see 
Section VIII]. • 

26. For example if is fully normalized so that 


/ <X (s) !x (n S ^ > “ <S n then: 

£ *q ’*o n,0 


b-0 

„( 0 ) _ 


.(I) 




rs 1 r 1 


27. 


C f ’ - - I < *f ) l^ U> - I a i + 1 a 2 • S "‘ 

<£ 3) ■ V 2 + l a 3 ! ' 

etc. Here the a^_ are constants which must be determined by some . ^ 
additional condition. Usually they are taken to be zero. 

t 
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